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ABSTRACT 

We  have  developed  a  new  method  for  determining  the  corotation  radii  of  density  waves  in  disk  galaxies,  which 
makes  use  of  the  calculated  radial  distribution  of  an  azimuthal  phase  shift  between  the  potential  and  density  wave 
patterns.  The  approach  originated  from  improved  theoretical  understanding  of  the  relation  between  the  morphology 
and  kinematics  of  galaxies  and  of  the  dynamical  interaction  between  density  waves  and  the  basic-state  disk  stars, 
which  results  in  the  secular  evolution  of  disk  galaxies.  In  this  paper  we  present  the  rationales  behind  the  method  and 
the  first  application  of  it  to  several  representative  barred  and  grand-design  spiral  galaxies,  using  near-infrared  images 
to  trace  the  mass  distributions,  as  well  as  to  calculate  the  potential  distributions  used  in  the  phase-shift  calculations. 
We  compare  our  results  with  those  from  other  existing  methods  for  locating  the  corotations  and  show  that  the  new 
method  both  confirms  the  previously  established  trends  of  bar-length  dependence  on  galaxy  morphological  types  and 
provides  new  insights  into  the  possible  extent  of  bars  in  disk  galaxies.  The  method  also  facilitates  the  estimation  of 
mass  accretion /excretion  rates  due  to  bar  and  spiral  density  waves,  providing  an  alternative  way  of  quantifying  the 
importance  of  these  features  in  disk  galaxies.  A  preliminary  analysis  of  a  larger  sample  shows  that  the  phase-shift 
method  is  likely  to  be  a  generally  applicable,  accurate,  and  essentially  model-independent  method  for  determining  the 
pattern  speeds  and  corotation  radii  of  single  or  nested  density  wave  patterns  in  galaxies.  Other  implications  of  the 
results  of  this  work  include  that  most  of  the  nearby  bright  disk  galaxies  appear  to  possess  quasi-stationary  spiral 
modes;  that  these  density  wave  modes,  as  well  as  the  associated  basic  states  of  the  galactic  disks,  slowly  transform 
over  the  time  span  of  a  Hubble  time  due  to  a  collective  dissipation  process  directly  related  to  the  presence  of  the  phase 
shift  between  the  potential  and  density  patterns;  and  that  self-consistent  A-particle  systems  contain  physics  not 
revealed  by  the  passive  orbit  analysis  approaches. 
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1.  INTRODUCTION 

More  than  two-thirds  of  all  bright  galaxies  in  the  nearby  uni¬ 
verse  are  disk  galaxies  (de  Vaucouleurs  1963;  Buta  et  al.  1994). 
The  fraction  of  disk-dominated  galactic  systems  also  increases 
steadily  toward  higher  redshift.  Understanding  the  kinematics, 
dynamics,  and  long-term  evolution  of  disk  galaxies  therefore  holds 
the  key  to  understanding  the  observed  morphological  transforma¬ 
tion  of  galaxies  across  the  Hubble  time  (Butcher  &  Oemler  1978a, 
1978b;  Dressier  1980;  Koo  &  Kron  1992;  Ellis  1997;  Couch  et  al. 
1994;  Dressier  et  al.  1994;  Lilly  et  al.  1998). 

The  most  prominent  feature  in  almost  all  disk  galaxies  is  the 
presence  of  global  patterns  such  as  bars,  rings,  and  spiral  struc¬ 
tures.  These  patterns  have  been  attributed  to  density  waves  which 
are  self-sustained  and  long-lasting  (Lin  &  Shu  1964).  Statistical 
evidence  indicates  that  most  of  these  patterns  can  survive  for  a 
significant  fraction  of  a  Hubble  time  (Elmegreen  &  Elmegreen 
1983,  1989);  A-body  simulations  of  spiral  disks  have  also  sub¬ 
stantiated  this  speculation  (Donner  &  Thomasson  1994;  Zhang 
1998).  These  long-lasting  density  wave  patterns  are  found  to  be, 
in  fact,  density  wave  modes ,  formed  by  the  feedback  interactions 
of  oppositely  propagating  density  wave  trains  in  the  radial  direc¬ 
tion,  with  an  amplification  mechanism  at  the  so-called  corotation 
resonance  (CR)  radius  where  the  density  wave  pattern  and  the  un¬ 
derlying  differentially  rotating  disk  matter  rotate  at  the  same  an¬ 
gular  speed  (Mark  1976;  Toomre  1981).  For  galaxies  which 


have  several  nested  spiral/bar  patterns,  it  is  expected  that  differ¬ 
ent  pattern  speeds  and  therefore  different  corotation  radii  exist 
for  the  different  patterns.  The  accurate  determination  of  the  pat¬ 
tern  speeds  and  corotation  radii  of  the  density  wave  patterns  is  a 
key  ingredient  in  determining  the  dynamical  state  of  the  parent 
galaxies,  which  will  help  address  their  evolutionary  history  and 
fate;  the  pattern  speeds  of  secondary  bar  patterns  also  provide 
clues  on  the  dark  matter  distributions  in  the  inner  region  and  thus 
have  cosmological  implications  (Debattista  &  Sellwood  2000). 

The  determination  of  the  spiral  density  wave  pattern  speed  and 
corotation  radius  has  historically  been  difficult,  despite  the  great 
variety  of  methods  proposed  in  the  past  to  measure  either  pa¬ 
rameter.  These  methods  are  based  on  various  physical  criteria  or 
assumptions.  The  earliest  proposal  for  corotation  determination 
was  based  on  Lin’s  (1970)  conjecture  that  an  organized  two¬ 
armed  spiral  pattern  exists  only  within  the  corotation  circle,  and 
arms  become  more  fragmented  outside.  Among  the  first  to  apply 
this  conjecture  to  physical  galaxies  were  Shu  et  al.  (1971),  who 
placed  the  corotation  radius  where  the  distribution  of  H  n  regions 
is  seen  to  end.  Tremaine  &  Weinberg  (1984,  hereafter  TW84)  de¬ 
veloped  an  approach  for  pattern  speed  determination  which  was 
based  on  the  assumption  of  continuity  of  a  luminosity  density 
tracer  and  the  measurement  of  surface  brightness  and  velocity 
distributions  of  the  tracer.  Canzian  (1993)  proposed  a  method  for 
corotation  determination  that  relied  on  the  qualitative  difference 
in  the  appearance  of  the  residual  velocity  field  of  a  density  tracer 
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inside  versus  outside  of  corotation.  Elmegreen  et  al.  (1989)  iden¬ 
tified  resonances  through  breaks  in  spiral  arms,  whereas  Elmegreen 
et  al.  (1998)  determined  corotation  from  the  transition  of  the  inward- 
to  outward-streaming  gas  velocity.  Buta  &  Purcell  (1998)  and  Buta 
et  al.  (1998)  used  the  locations  of  outer  rings  to  infer  the  location  of 
corotation  in  several  resonance  ring  galaxies.  Clemens  &  Alexander 
(2001)  used  high-resolution  spectroscopic  data  to  find  signatures 
of  the  preshock  and  postshock  gas  across  the  spiral  arms  and 
used  this  information  to  determine  the  shock  velocity  and  hence 
the  pattern  speed.  Oey  et  al.  (2003)  used  an  H  n  region  isochrone 
technique  to  constrain  the  corotation  radius.  Various  numerical 
simulation  techniques  (Hunter  et  al.  1988;  England  1989;  Ball 
1992;  Lindblad  et  al.  1996;  Lindblad  &  Kristen  1996;  Salo  et  al. 
1999;  Rautiainen  et  al.  2005)  attempted  to  visually  match  ob¬ 
served  galaxy  morphology  with  the  simulated  density  wave  of  a 
given  pattern  speed,  using  mostly  non-self-consistent  (i.e.,  en¬ 
forced)  potential  waves.  Other  methods  used  offsets  between  Ha 
and  CO  arms  (Egusa  et  al.  2004)  or  between  B-  and  /-band  spi¬ 
ral  arms  (Puerari  &  Dottori  1997)  to  determine  the  corotation 
radii.  Almost  all  of  the  above-mentioned  methods  have  only 
limited  ranges  of  applicability;  many  lack  either  accuracy  or  model 
independency. 

In  what  follows  we  describe  a  novel  and  general  approach  to 
locating  the  corotation  radii  of  density  waves  in  disk  galaxies, 
based  on  new  theoretical  understandings  of  the  properties  of  quasi¬ 
stationary  density  wave  modes  (Zhang  1996,  1998,  1999).  The 
technique  employs  the  use  of  near-infrared  (NIR)  images  from 
ground-  and  space-based  surveys  to  infer  densities  and  gravi¬ 
tational  potentials.  The  NIR  region  of  the  spectrum  contains  ra¬ 
diation  contributions  mainly  from  old  giant  and  supergiant  stars 
(Frogel  et  al.  1996)  and  is  considered  a  much  better  tracer  of  the 
stellar  mass  distribution  than  optical  bands.  In  the  NIR  the  effects 
of  dust  extinction  are  also  considerably  reduced. 

The  application  of  new  theoretical  results  of  the  relation  be¬ 
tween  the  morphological  characteristics  and  the  kinematic  states 
of  galaxies  provides  both  a  practical  approach  for  corotation  deter¬ 
mination  and  a  direct  quantitative  comparison  between  theoretical 
predictions  and  observed  galaxy  properties.  A  further  applica¬ 
tion  of  the  same  general  approach  is  the  estimation  of  the  mass 
accretion /excretion  rates  due  to  bars  and  spirals.  Apart  from 
providing  a  new  way  of  quantifying  the  importance  of  these  fea¬ 
tures  in  galaxies  (see  Buta  et  al.  2005  and  references  therein  for 
other  approaches),  knowledge  of  these  rates  is  important  to  our 
understanding  of  the  secular  evolutionary  transformation  of  gal¬ 
axy  morphological  types  (Zhang  &  Buta  2007). 


sity  wave  patterns.  These  global  patterns  were  shown  to  drive  a 
secular  energy  and  angular  momentum  exchange  process  between 
the  disk  matter  and  the  wave  pattern,  mediated  by  a  local  gravi¬ 
tational  instability,  or  a  collisionless  shock,  at  the  potential  mini¬ 
mum  of  the  pattern  (Zhang  1996).  The  integral  manifestation  of 
this  process  is  an  azimuthal  phase  shift  between  the  potential  and 
density  spirals,  which  results  in  a  secular  torque  action  between 
the  wave  pattern  and  the  underlying  disk  matter.  As  a  result  of  the 
torqueing  of  the  wave  on  the  disk  matter,  the  matter  inside  the  CR 
radius  loses  angular  momentum  to  the  wave  secularly  and  sinks 
inward.  The  wave  carries  the  angular  momentum  it  receives  from 
the  inner  disk  matter  to  the  outer  disk  and  deposits  it  there,  causing 
the  matter  in  the  outer  disk  to  drift  further  out.  The  resulting 
secular  morphological  evolution  causes  the  Hubble  type  of  an 
average  galaxy  to  evolve  from  late  to  early  (Zhang  1999). 

One  of  the  predictions  of  the  above  theoretical  work  is  that  for 
a  self-sustained  spiral  or  bar  mode,  the  radial  distribution  of  the 
potential-density  phase  shift  should  be  such  that  it  changes  sign 
at  the  corotation  radius  (Zhang  1996, 1998).  (Note  that  the  phase- 
shift  expression  in  eq.  [54]  of  Zhang  [1996]  formally  diverges 
exactly  at  corotation  due  to  the  breakdown  of  the  orbit  approach 
at  the  exact  resonance,  but  Zhang  [1998]  shows,  through  the  re¬ 
lation  of  phase  shift  and  angular  momentum  flux  distribution,  that 
for  a  spontaneously  formed  density  wave  mode  the  phase  shift 
crosses  zero  exactly  at  corotation.)  This  happens  because  the  wave 
rotates  more  slowly  than  the  disk  stars  inside  corotation  (and  thus 
the  wave  has  negative  energy  and  angular  momentum  density 
with  respect  to  the  disk  stars)  and  vice  versa  outside  corotation.  A 
positive  sign  of  the  phase  shift,  defined  as  when  the  potential 
pattern  lags  the  density  pattern  in  azimuth,  for  a  radial  location 
inside  r(CR),  leads  on  average  to  an  angular  momentum  loss  from 
the  population  of  disk  stars  in  a  given  annulus  inside  corotation 
(since  the  phase  shift  was  shown  to  be  directly  correlated  with  the 
gradient  of  the  angular  momentum  flux  being  carried  outward  by 
the  wave),  which  further  leads  to  the  spontaneous  growth  of  the 
wave  mode  in  the  linear  regime  and  to  the  damping  of  the  growing 
wave  amplitude  toward  its  quasi- steady  state  value  in  the  non¬ 
linear  regime  (Zhang  1998). 

In  the  quasi- steady  state,  the  rate  of  angular  momentum  ex¬ 
change  between  a  skewed  density  wave  pattern  and  the  basic  state1 
of  the  disk  per  unit  area,  due  to  the  torquing  of  the  wave  on  the  disk, 
is  given  by  (Zhang  1996,  1998) 
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2.  POTENTIAL-DENSITY  PHASE-SHIFT  METHOD 
FOR  COROTATION  DETERMINATION 

Our  quantitative  investigation  of  the  kinematic  and  dynamical 
properties  of  galaxies  using  NIR  images  is  motivated  by  new 
theoretical  insight  into  the  role  of  density  waves  in  the  long-term 
maintenance  and  evolution  of  galaxies.  Although  density  wave 
theory  was  developed  more  than  40  years  ago  (Lin  &  Shu  1964), 
there  has  been  for  most  of  this  time  a  commonly  held  belief  that  the 
stellar  motion  in  a  galaxy  containing  quasi-stationary  wave  patterns 
conserves  the  Jacobi  integral  in  the  rotating  frame  of  the  pattern 
(Binney  &  Tremaine  1987,  eq.  [3-88]).  The  orbital  radii  of  stars 
under  such  applied  potentials  generally  will  not  exhibit  secular 
decay  or  increase,  and  there  is  no  wave  or  basic-state  (i.e.,  the  axi- 
symmetric  part  of  the  disk)  interaction  except  at  the  wave-particle 
resonances  (Lynden-Bell  &  Kalnajs  1972,  hereafter  LBK72). 

It  was  first  demonstrated  in  Zhang  (1996,  1998,  1999)  that 
secular  orbital  changes  of  stars  across  the  entire  galaxy  disk  are 
in  fact  possible  due  to  a  collective  instability  induced  by  the  den¬ 


where  £i  represents  the  perturbation  density  waveform  and  Vi  is 
the  perturbation  potential  waveform  (with  the  circularly  sym¬ 
metric  m  =  0  component  already  subtracted).  For  two  sinusoidal 
waveforms,  equation  (1)  can  be  further  written  as 

dL  f  m  \ 

-pr)  =  (—  JAs(r)Av(r)  sin  [m<j>o(r)],  (2) 

where  and^y  are  the  amplitudes  of  the  density  and  potential 
waves,  respectively,  m  is  the  number  of  spiral  arms,  and  is 
the  phase  shift  between  these  two  waveforms,  defined  as  being 


1  The  term  “basic  state”  has  been  used  to  describe  the  properties  of  the  axi- 
symmetric  part  of  the  disk  which  forms  the  boundary  conditions  for  the  pertur¬ 
bative  studies  of  the  nonaxisymmetric  density  waves  and  modes  (see,  e.g.,  Lin  & 
Lau  1979).  These  properties  usually  include  the  radial  distribution  of  surface  den¬ 
sity  of  the  disk,  the  radial  distribution  of  the  circular  velocity  (which  incorporates 
the  influence  of  any  rigid  bulge  and  halo  component),  and  the  radial  distribution 
of  velocity  dispersion. 
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positive  if  the  potential  lags  the  density  in  the  azimuthal  direction 
in  the  sense  of  the  galactic  rotation.  For  most  galaxies  the  sense 
of  galactic  rotation  can  be  determined  by  the  assumption  that  the 
density  wave  pattern  is  trailing,  although  exceptions  do  exist 
which  will  need  to  be  studied  on  a  case-by-case  basis. 

Using  equations  (1)  and  (2)  we  see  that  the  radial  distribution 
of  an  equivalent  azimuthal  phase  shift  </>o(r,  <j>)  between  the  (gen¬ 
erally  nonlinear  and  nonsinusoidal)  potential  and  density  pat¬ 
terns  for  spirals  or  bars  in  a  disk  galaxy  can  be  calculated  from 


0o(r,  (/))  =  —  sin 
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1  jfyi(gVi/g^ 

m  9  ’ 


where  the  potential  V  is  calculated  from  the  Poisson  integral  of 
the  density  distribution,  with  the  density  distribution  itself  calcu¬ 
lated  from  the  NIR  image  of  a  galaxy  for  an  assumed  mass-to-light 
ratio  {MIL)  distribution.  The  equivalent  phase  shift  is  the  amount 
of  phase  shift  which  would  be  present  between  two  sinusoidal 
waveforms  if  each  were  endowed  with  the  same  energy  as  the 
corresponding  nonlinear  waveform  and  which  would  lead  to  the 
same  value  for  the  torque  integral  as  would  the  nonlinear  wave¬ 
forms.  The  phase  shift  is  in  general  nonzero  as  long  as  the  density 
pattern  is  skewed,  i.e.,  in  the  case  of  spirals,  twisted  or  offset  bars, 
or  even  some  twisted  three-dimensional  mass  distributions  as  ob¬ 
served  in  many  high-redshift  protogalaxies. 

Note  that  even  though  these  waveforms  are  given  in  their 
perturbational  form  (with  subscript  1  in  the  above  equations),  in 
carrying  out  the  actual  calculations  we  can  simply  use  the  non- 
perturbed  quantities  in  the  numerator,  since  the  effect  of  the  cir¬ 
cularly  symmetric  components  in  V  and  E  will  differentiate  and 
integrate  out  to  zero.  This  reduces  the  sensitivity  to  uncertainties 
in  the  MIL  of  any  axisymmetric  mass  component  (or,  in  fact,  any 
nonskewed  component,  including  triaxial  luminous  or  dark  halo 
components),  including  the  axisymmetric  bulge,  halo,  and  disk 
components.  The  normalization  factor  in  the  denominator,  how¬ 
ever,  will  be  affected  by  the  axisymmetric  components,  so  these 
should  be  subtracted  out  before  doing  the  normalization.  However, 
any  uncertainty  in  the  normalization  has  no  effect  on  the  locations 
of  the  positive-to-negative  crossings  of  the  phase  shift  versus  ra¬ 
dius  curve,  which  are  what  we  use  in  the  corotation  determination. 

A  practical  note  for  the  sign  convention  of  the  phase  shift  is 
that  we  assume  the  phase  shift  to  be  positive  when  the  potential 
pattern  lags  the  density  pattern  in  the  direction  of  galactic  rota¬ 
tion.  Since  the  face-on  view  of  the  spiral  (or  twisted  bar)  disk  comes 
in  two  flavors,  i.e.,  either  the  S-shaped  pattern  or  the  Z-shaped 
pattern,  these  in  general  indicate  a  counterclockwise  or  clock¬ 
wise  rotation  direction,  respectively,  if  the  pattern  is  assumed  to 
be  trailing.  In  the  above  definitions,  therefore,  the  sense  of  the 
phase  shift  is  correct  only  for  the  S-shaped  spiral  or  bar  (which 
winds  and  rotates  in  the  same  counterclockwise  direction  as  used 
in  the  azimuthal  angle  definition).  For  a  Z-shaped  pattern  the  cal¬ 
culated  phase-shift  curve  should  be  negated  to  obtain  the  correct 
sense  of  the  phase  shift  in  the  direction  of  galactic  rotation.  In 
practice,  we  usually  generate  both  versions  (called  a  “plus”  plot 
and  a  “minus”  plot)  for  the  phase  shift  versus  radius  distributions 
and  compare  both  with  the  image  of  the  galaxy,  since  in  certain 
cases  the  sense  of  the  winding  is  not  so  obvious  (for  some  barred 
galaxies  where  the  bars  are  fairly  straight  or  for  cases  where 
leading  arms  may  be  present,  as  for  the  case  of  NGC  4622,  which 
we  discuss  in  §  3.4). 

After  obtaining  the  phase  shift  versus  galactic  radius  plot,  the 
positions  of  the  successive  corotation  radii  can  be  read  off  as  the 


positive-to-negative  crossings  of  the  phase-shift  plot.  These  mark 
the  locations  where  the  direction  of  angular  momentum  flow  be¬ 
tween  the  disk  matter  and  the  density  wave  changes  sign,  at  the 
quasi-steady  state  of  the  wave.  The  negative-to -positive  cross¬ 
ings  of  the  phase  shift,  on  the  other  hand,  generally  mark  the 
transition  locations  of  nested  modal  structures;  i.e.,  each  such 
negative-to-positive  crossing  is  likely  to  be  where  the  pattern 
speeds  of  the  density  waves  change  abruptly.  This  is  because  the 
density  wave  should  have  positive  angular  momentum  inside  co¬ 
rotation  and  negative  angular  momentum  outside  corotation.  The 
negative-to-positive  crossings  mark  the  boundaries  where  the 
wave  deposits  angular  momentum  onto  the  disk  (for  the  inner  mode 
outside  its  corotation)  or  the  wave  takes  away  angular  mo¬ 
mentum  (for  the  outer  mode  inside  its  own  corotation).  So  the 
phase-shift  distribution  provides  not  only  a  way  for  the  objective 
determination  of  density  wave  corotation  radii,  but  also  a  way 
for  the  objective  determination  of  the  radial  extent  of  the  indi¬ 
vidual  modes  in  the  nested-mode  cases;  i.e.,  the  effective  radii 
where  adjacent  modes  decouple,  something  that  methods  such  as 
TW84  could  not  determine  in  a  model-independent  way.  This 
assertion,  of  course,  is  based  on  the  assumption  that  the  radial 
penetration  of  the  different  sets  of  modes  are  minimal,  which 
should  be  expected  if  the  morphology  of  the  modes  is  quasi¬ 
steady,  since  the  interpenetration  of  modes  with  different  pat¬ 
tern  speeds  will  result  in  rapidly  changing  overall  density  wave 
morphology. 

We  emphasize  that  for  the  above  approach  to  work,  the  indi¬ 
vidual  density  wave  modes  in  the  galaxies  of  interest  have  to  have 
already  achieved  a  quasi-stationary  state ,  since  otherwise  one 
would  not  expect  the  purely  morphological  characteristics  (i.e., 
the  image  of  a  galaxy)  to  correlate  with  kinematic  status  informa¬ 
tion  (i.e.,  the  rotational  state  of  the  disk  stars  and  density  waves). 
Therefore,  a  demonstration  of  the  validity  of  this  approach,  i.e., 
whether  the  corotations  we  determined  actually  correspond  to 
visible  resonant  features  in  the  images  of  galaxies,  also  confirms 
the  quasi- steady  state  of  the  density  wave  modes  in  the  relevant 
galaxies. 

3.  APPLICATION  TO  INDIVIDUAL  GALAXIES 

The  potential-density  phase  shift,  as  well  as  its  sign  change 
across  corotation,  has  been  confirmed  with  TV-body  simulations 
(Zhang  1996, 1998),  which  verified  that  the  method  gives  the  same 
corotation  radius  as  that  derived  from  the  pattern  speed  of  the  den¬ 
sity  wave  in  conjunction  with  a  galactic  rotation  curve.  Our  present 
study  represents  the  first  detection  of  the  phase  shift  in  observed 
galaxies,  based  mainly  on  //-band  (1.65  gm)  images  from  the 
Ohio  State  University  Bright  Galaxy  Survey  (OSUBGS;  Eskridge 
et  al.  2002),  as  well  as  on  other,  smaller  samples. 

NIR  images  can  be  used  to  measure  the  phase  shifts  because, 
as  previously  noted,  such  images  trace  the  stellar  mass  distri¬ 
bution  better  than  optical  images,  and  for  this  reason  they  are 
also  the  best  to  use  for  calculating  gravitational  potentials  (e.g., 
Quillen  et  al.  1994).  For  most  of  the  calculations  given  below, 
the  MIL  has  been  assumed  to  be  a  constant  independent  of  radius 
unless  otherwise  noted.  The  exact  value  of  this  constant  is  not 
important,  since  it  cancels  out  between  the  numerator  and  the 
denominator  in  the  definition  of  equation  (3).  The  validity  of  this 
assumption  is  further  justified  in  §  3.2.  Following  Quillen  et  al., 
we  calculate  the  two-dimensional  (2D)  potential  at  the  plane  of 
the  galaxy  from 


$(x,y,z)  =  -G 


j  H(x'  ,y')g(x  —  x'  ,y  —  y')dx'  dy' , 


(4) 
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where  the  2D  Green’s  function  g(r)  =  g{\/x2  +  y2)  is  given  by 


p(z)  dz 
yjr2  +  z2  ’ 


(5) 


where  p(z)  is  the  normalized  z-distribution  of  the  volume  density 
of  matter,  assumed  to  be  independent  of  galactic  radii  r.  In  prac¬ 
tice,  we  find  that  the  several  different  forms  of  the  assumed  p(z) 
distribution  given  in  Quillen  et  al.  (1994)  do  not  give  significantly 
different  phase-shift  results,  especially  for  the  zero  crossings,  even 
though  they  do  change  the  absolute  values  of  the  calculated  po¬ 
tential  by  about  20%.  This  insensitivity  to  the  exact  z-profile  of  the 
density  distribution  is  related  to  the  fact  that  phase  shift  is  deter¬ 
mined  mostly  by  the  global  distribution  of  the  pattern  pitch  angle 
and  radial  density  variation  and  is  insensitive  to  how  “puffed  up” 
the  pattern  is  in  the  z-direction.  We  have  available  Quillen  et  al.  ’s 
original  code  for  potential  calculation  written  in  C.  But  since  that 
code  assumes  arbitrary  units,  which  was  fine  for  the  phase  shift 
calculation  but  not  adequate  for  the  mass-accretion  calculation  in 
which  we  are  also  interested  (see  §  4.3),  we  completely  rewrote 
the  code  in  FORTRAN  and  adopted  a  different  fast  Fourier  trans¬ 
form  routine.  The  resulting  code  gives  potential  in  physical  units 
and  consistently  outperforms  Quillen  et  al.  ’s  original  code  in 
terms  of  resolving  the  central  structures  of  galaxies  in  the  phase- 
shift  results  (when  the  imaging  data  used  have  high  enough 
spatial  resolution),  which  turned  out  to  be  important  for  galaxies 
possessing  nested  density  wave  modes,  such  as  NGC  1530  and 
NGC  4321.  An  example  of  the  calculated  potential  for  NGC 
1530  using  the  above  approach  is  given  in  §  3.2. 

After  the  surface  density  map  and  the  potential  map  for  each 
galaxy  are  thus  determined  (on  the  Cartesian  grid),  the  phase  shift 
versus  galactic  radius  curve  is  calculated  using  equation  (3)  on  an 
exponential  polar  grid  similar  to  that  used  in  TV-body  simula¬ 
tions  (Zhang  1 996),  with  values  of  the  potential  and  density  on  the 
polar  grid  obtained  through  an  interpolation  between  the  polar  and 
Cartesian  grid.  The  azimuthal  distributions  for  the  potential  and 
density  at  each  radius  on  the  polar  grid  are  first  fitted  with  cubic 
spline  curves;  the  differentiation  and  integration  are  then  performed 
based  on  the  fitted  potential  and  density  curves  using  routines  in 
Press  et  al.  (1992).  For  most  galaxies  the  convergence  of  the  inte¬ 
gration  is  rapid.  Occasionally,  however,  when  the  phase  shift  versus 
radius  curve  is  particularly  noisy,  there  exist  some  radii  where  the 
phase-shift  calculations  fail  to  converge  after  a  reasonable  number 
of  iterations  (so  far  this  happens  for  less  than  1%  of  the  galaxies  we 
have  calculated).  Since  this  occurs  for  isolated  radii  only  we  used 
interpolation  of  the  phase-shift  results  of  neighboring  points  to 
arrive  at  the  values  for  these  radial  locations. 

In  the  rest  of  this  section,  we  describe  representative  results  of 
the  phase-shift  calculation  and  corotation  determination  for  a 
number  of  galaxies  we  have  analyzed  so  far,  and  we  defer  the 
description  of  such  calculations  for  the  majority  of  the  galaxies 
in  the  OSUBGS  sample  to  a  sequel  paper  (R.  J.  Buta  &  X.  Zhang 
2007,  in  preparation).  The  detailed  assessment  of  the  accuracy 
of  the  method,  including  the  dependence  of  the  results  on  the 
uncertainties  in  vertical  scale  height,  the  extent  of  the  galaxy 
image  size  used,  the  signal-to-noise  ratio  of  the  image,  and  the 
wave  band  of  the  image  used,  as  well  as  on  the  uncertainties  in 
the  orientation  parameters  used  for  the  deprojection  of  the  im¬ 
age,  will  also  be  given  there. 

3.1.  Representative  Cases 

Among  the  galaxies  we  have  analyzed  so  far,  the  simplest 
cases  to  interpret  (and  often  also  the  most  striking  in  appearance) 
are  those  where  the  phase  shift  versus  galactic  radius  curve  has 


only  a  single  major  positive-to-negative  crossing  across  the 
entire  disk,  although  some  may  contain  a  second,  or  even  a  third, 
minor  crossing  in  the  central  region  of  the  galaxy  as  well,  sig¬ 
naling  secondary,  nested  nuclear  spiral  or  bar  patterns.  Examples 
from  the  OSUBGS  which  fall  into  this  class  include  NGC  4051, 
4314,  4548,  4665,  5054,  5247,  5248,  6215,  7479,  and  7552. 

The  galaxies  listed  above  can  be  further  divided  into  two  sub¬ 
categories:  (1)  those  with  prominent  bars,  such  as  NGC  4314, 
4548,  and  7479,  which  have  a  derived  CR  circle  that  lies  very  near 
to  (or  slightly  farther  out  than)  the  end  of  the  bar;  and  (2)  those 
containing  grand-design  spirals  such  as  NGC  5054,  5247,  and 
5248,  where  a  single  CR  circle  lies  in  the  middle  (or  sometimes 
near  the  end)  of  the  main  spiral  pattern. 

For  the  kind  of  bars  where  r(CR)  is  at  the  end  of  the  bar  and 
there  is  a  spiral  pattern  emanating  from  the  bar  end,  the  inner  and 
outer  patterns  are  most  likely  corotating;  i.e.,  in  the  case  in  which 
the  bar  ends  at  its  r(CR)  we  most  likely  have  a  bar-driven  spiral 
outside  of  the  bar.  For  the  case  of  a  spiral  pattern  with  a  single 
r  (CR)  ending  in  the  middle  of  the  spiral  pattern,  the  spiral  is  likely 
to  have  extended  all  the  way  to  its  outer  Lindblad  resonance 
(OLR). 

In  Figure  1  (left)  we  show  the  calculated  phase  shift  versus  ga¬ 
lactic  radius  plot  for  the  nearly  face-on  SBa  galaxy  NGC  4314, 
based  on  a  deprojected  OSUBGS  H- band  image.  The  gravita¬ 
tional  potential  has  been  calculated  under  the  assumption  of  a 
constant  MIL  and  an  exponential  vertical  density  distribution 
with  scale  height  hz  of  one-fourth  the  radial  scale  length  Hr  (de 
Grijs  1998;  Laurikainen  et  al.  2004),  using  the  approach  described 
before.  The  location  of  the  corotations  (two  in  this  case)  are  marked 
by  the  positive-to-negative  crossings  of  the  phase-shift  curve  (ar¬ 
rows)  at  9.0"  and  76.5".  In  Figure  1  (right)  these  two  corotation 
radii  are  plotted  on  the  same  H- band  image. 

As  can  be  seen  from  the  image  for  NGC  4314  overlaid  with 
the  corotation  circles,  the  outer  corotation  corresponds  to  the  end 
of  the  bar.  This  outer  corotation  radius  determined  by  the  phase- 
shift  method  is  found  to  be  very  close  to  the  value  determined  for 
this  galaxy  by  Quillen  et  al.  (1994).  The  inner  corotation  indicates 
the  existence  of  an  inner  density  wave  pattern,  likely  associated 
with  the  well-known  nuclear  ring/spiral  in  this  galaxy  (Sandage 
1961).  This  feature  has  a  radius  of  8",  according  to  Buta  & 
Crocker  (1993).  Note  that  even  though  the  inner  density  wave 
pattern  is  not  well  resolved  due  to  the  limited  spatial  resolution 
of  this  H- band  image,  a  potential-density  phase-shift  crossover  is 
nonetheless  clearly  predicted.  This  apparent  “superresolution” 
effect  is  achieved  because  the  potential  is  obtained  as  the  Poisson 
integral  of  the  global  density  distribution  of  the  entire  disk  matter 
and  is  not  determined  solely  by  the  local  density  distribution  near 
the  central  region  of  the  galaxy. 

In  Figure  2  we  show  the  classic  case  of  an  ordinary  grand- 
design  spiral  pattern,  NGC  5247,  and  the  corotation  circles  deter¬ 
mined  through  the  phase-shift  method.  The  outer  crossing,  which 
indicates  the  location  of  the  main  CR  circle,  is  in  the  middle  of  the 
outer  spiral  arms.  The  outer  corotation  radius  we  determined  ap¬ 
pears  to  be  close  to  where  the  inner  dust  lanes  tmncate  on  the  two 
spiral  arms  (this  effect  is  more  obvious  from  a  5-band  image  of 
this  galaxy),  which  is  usually  a  good  and  reliable  indication  of  the 
corotation  radius. 

Contopoulos  &  Grosbol  (1986,  1988)  had  previously  deter¬ 
mined  the  corotation  radius  for  this  galaxy  under  the  assumption 
that  the  4  : 1  resonance  (where  D  —  ft/4  =  Llp,  where  ft  is  the  epi¬ 
cycle  frequency,  D  is  the  angular  speed  of  the  disk  matter,  and  Ctp 
is  the  pattern  speed  of  the  density  wave)  lies  near  the  end  of  the 
strong  spiral  pattern.  Their  analysis  placed  this  resonance  at  12  kpc, 
which  corresponds  to  116"  assuming  a  distance  of  21.3  Mpc. 
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Fig.  1 . — Left.  Calculated  phase  shift  vs.  radius  for  OSUBGS  galaxy  NGC  43 14,  type  SBa.  The  arrows  point  to  the  locations  of  two  positive-to-negative  phase-shift  crossings 
(radii  9.0"  and  76.5"),  interpreted  as  CR  radii  of  the  inner  and  outer  patterns.  The  short  horizontal  line  indicates  the  radial  range  of  the  bar  obtained  based  on  two  definitions 
given  in  §  3.5.  Right:  Deprojected  NIR  //-band  (1 .65  gm)  image  of  NGC  43 14  in  log  scale,  with  the  corotations  determined  by  the  phase-shift  method  superimposed  {circles). 


The  corresponding  location  of  the  corotation  radius  from  their 
analysis  is  23  kpc,  or  222",  well  outside  the  NIR  image  of  the 
spiral  pattern.  We  note  that  Contopoulos  &  Grosbol’s  original 
analysis  was  based  on  matching  the  surface  density  obtained  from 
passive  orbital  responses  calculated  under  a  forced  potential  to 
that  of  the  given  galaxy  morphology.  They  also  assumed  that  the 
sharp  truncation  of  the  spiral  pattern  occurs  near  the  4  : 1  reso¬ 
nance,  as  the  starting  point  of  their  model  construction  process. 
We  present  extensive  arguments  (in  §  4.1)  of  why  the  passive 
orbit  approach  for  studying  the  extent  of  density  waves  in  gal¬ 
axies  has  limitations  when  collective  dissipation  sets  in,  as  is  the 
case  for  the  density  wave  modes  in  physical  galaxies. 

3.2.  Nested  Resonances  and  Effect  of  M/L 

NGC  1530  is  a  high-luminosity  SB(rs)b  spiral  with  one  of  the 
strongest  known  bars.  We  use  a  Ks -band  (2.15  /mi)  image  of  this 
galaxy  (Block  et  al.  2004)  to  conduct  more  detailed  analysis  of 
the  nested  resonance  features  already  encountered  for  the  two 
galaxies  in  §  3.1.  We  also  use  this  galaxy  as  an  example  to 


illustrate  the  effect  of  MIL  (and  its  uncertainties)  on  the  phase 
shift  and  corotation  determination,  as  well  as  on  determining  the 
mass  flow  rate  (accretion/excretion)  as  a  function  of  radius  (§  4.3). 
A  detailed  image  of  the  calculated  potential  for  this  galaxy  is 
given  as  a  representative  example,  using  the  numerical  codes  we 
described  at  the  beginning  of  §  3. 

We  examine  MIL  variations  using  the  color-dependent  for¬ 
mulae  of  Bell  &  de  Jong  (2001).  The  conversion  of  the  i^-band 
image  of  NGC  1530  into  a  surface  mass  density  map  was  accom¬ 
plished  using  a  V  —  Ks  color  index  map.  The  f^band  image  was 
obtained  with  the  2.5  m  Nordic  Optical  Telescope,  while  the  Ks- 
band  image  was  obtained  with  the  4.2  m  William  Herschel  Tele¬ 
scope.  These  images  were  matched  in  scale  and  orientation 
using  IRAF  routines  GEOMAP  and  GEOTRAN.  Zero  points 
were  obtained  using  published  photoelectric  aperture  photome¬ 
try  for  the  Hband  image  and  the  2MASS  A^-band  magnitude  in  a 
14"  aperture  provided  on  the  NED  Web  site.  Each  image  was 
cleaned  of  foreground  stars,  deprojected  using  a  position  angle 
of  8°  and  inclination  of  45°  from  Regan  et  al.  (1996),  and  rotated 


Fig.  2. — Left :  Phase  shift  vs.  radius  for  the  ordinary  spiral  NGC  5247,  showing  two  major  positive-to-negative  crossings  at  r  =  87"  and  14".  Right :  Deprojected  NIR 
//-band  (1.65  gm)  image  of  NGC  5247  in  log  scale,  with  the  corotations  determined  by  the  phase-shift  method  superimposed  {circles). 
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Fig.  3. — Left:  Radial  dependence  of  the  azimuthally  averaged  MIL  for  NGC  1530.  Right :  Interpolated  MIL  profile  along  the  bar  ( solid  curve )  and  90°  from  the  bar 
( dashed  curve). 


so  that  the  bar  was  horizontal.  The  resulting  V  —  Ks  color  index 
map  showed  strong  red  dust  lanes  in  the  bar  and  bluer  spiral  arms. 
Assuming  a  reference  bar  color  of  (V  —  ATs.)rcf  =  3.47  (corrected 
for  Galactic  extinction  as  given  by  NED  [Schlegel  et  al.  1998]), 
we  were  able  to  effectively  correct  the  color  index  map  for  the  dust 
lanes  using  a  simple  screen  approximation.2 

The  corrected  color  index  map  was  then  used  to  derive  a  MIL 
map  using  log  (MIL)  =  —1.087  +  0.314(E  —  K )  from  Table  1 
of  Bell  &  de  Jong  (2001),  assuming  V  —  K  «  V  —  Ks.  This  is 
for  the  “formation  epoch  model  with  bursts”  that  uses  a  scaled 
Salpeter  initial  mass  function,  and  it  gives  MIL  =  1.0  for  the 
reference  color  of  3.47.  Using  the  NED  zero  point,  the  X^-band 


2  The  reference  bar  color  refers  to  the  reddening-free  color.  Any  bar  point  redder 
than  this  is  assumed  to  be  reddened.  Any  galaxy  point  bluer  is  assumed  to  have  that 
color  as  intrinsic.  This  may  not  be  correct  in  every  case,  but  no  other  large-scale 
reddening  in  NGC  1530  is  as  serious  as  that  in  the  bar  dust  lanes.  The  screen 
approximation  assumes  that  all  the  dust  doing  the  reddening  is  in  the  foreground  of 
the  stars  being  reddened,  rather  than  mixed  with  the  stars.  Thus,  the  dust  acts  like  a 
“screen.”  This  is  not  perfectly  accurate,  but  it  worked  well  enough  for  our  analysis. 


image  was  converted  to  solar  X-luminosities  per  square  parsec 
using  an  absolute  magnitude  Mk(  0 )  =  3.33  from  Worthey 
(1994).  Multiplying  the  original  image  by  the  MIL  map  gave  the 
surface  mass  density  image.  The  main  effect  of  the  MIL  correc¬ 
tion  was  to  reduce  the  surface  density  of  the  spiral  compared  to  the 
bar,  thereby  weakening  the  effect  of  the  spiral.  The  azimuthally 
averaged  MIL  for  NGC  1530  is  given  in  Figure  3  (left),  along 
with  the  breakdown  of  the  MIL  along  the  bar  and  90°  from  the 
bar  (Fig.  3,  right).  This  procedure  gave  the  surface  mass  density 
map  in  Figure  4  (right),  which  is  in  units  of  M0  pc  “2. 

Figure  4  (left)  shows  the  calculated  phase  shift  versus  radius 
for  NGC  1530,  using  either  a  constant  MIL  (dashed  curve)  or  a 
position-dependent  MIL  (solid  curve).  We  observe  that  the  effect 
of  varying  the  MIL  changed  mostly  the  scalings  of  the  phase- 
shift  results,  and  the  zero  crossings  changed  by  only  about  1%.  This 
is  to  be  expected  from  the  arguments  given  in  §  2,  which  showed  that 
the  radial  MIL  change  for  the  axisymmetric  components  does  not 
impact  the  values  of  the  zero  crossings  of  the  phase  shift,  and 
much  of  the  MIL  variation  in  physical  galaxies  is  due  to  the  MIL 


r(kpc) 


Fig.  4. — Left :  Phase  shift  vs.  galactic  radius  calculated  for  the  SBb  galaxy  NGC  1530  from  a  AT*. -band  image.  The  solid  curve  was  calculated  using  a  position-dependent 
MIL  map,  and  the  dashed  curve  was  calculated  using  the  surface  density  map  assuming  constant  MIL.  Two  CR  circles  are  indicated.  The  short  horizontal  line  indicates  the 
range  for  the  end  of  the  bar  as  determined  based  on  the  two  definitions  given  in  §  3.5.  Right :  Deprojected  AC -band  image  (Block  et  al.  2004)  of  NGC  1530  in  log  scale,  with 
superimposed  corotation  (circles),  rotated  such  that  the  bar  is  horizontal.  This  image  has  been  corrected  for  stellar  MIL  variations  following  Bell  &  de  Jong  (2001).  The 
radius  scale  is  based  on  a  distance  of  36.6  Mpc  from  Tully  (1988). 
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change  of  the  axisymmetric  components  only,  since  the  popula¬ 
tion  of  stars  supporting  the  density  wave  is  relatively  homoge¬ 
neous  in  composition  (i.e.,  mostly  old  giant  and  supergiant  stars). 
Figure  4  {right)  displays  the  converted  surface  density  map  of 
the  galaxy  image  in  log  scale,  with  two  superimposed  corotation 
circles  determined  using  the  phase-shift  method.  The  radial  scalings 
used  for  NGC  1530  in  the  above  calculations  are  based  on  a  dis¬ 
tance  of  36.6  Mpc  from  Tully  (1988). 

For  reference,  we  also  present  in  Figure  5  an  image  of  the  cal¬ 
culated  potential  distribution  used  in  the  phase-shift  calculations 
for  this  galaxy.  From  this  figure  one  can  discern  the  slow  iso- 
photal  twist,  with  a  level  of  skewness  that  is  more  gradual  than 
the  density  pattern  skewness  observed  in  Figure  4  {right),  consis¬ 
tent  with  the  sense  of  “phase  shift  sign  change”  inside  and  out¬ 
side  corotation  (which  would  mean  that  a  potential  pattern  is  in 
general  straighter  than  the  density  pattern). 

As  was  the  case  for  NGC  43 14,  the  inner  CR  for  NGC  1 530  is  as¬ 
sociated  with  a  small  central  pattern  with  a  size  of  about  1 . 1  kpc, 


and  the  overall  extent  of  this  central  feature  appears  to  be  delineated 
by  the  negative-to-positive  crossing  on  the  phase-shift  plot.  Im¬ 
mediately  following  the  central  pattern  is  the  intermediate  oval 
pattern  between  1  and  2  kpc,  with  its  own  CR  at  about  1.3  kpc 
(this  was  marked  as  “CRi”  in  this  plot,  since  the  innermost  co¬ 
rotation  radius  is  not  very  well  resolved  by  the  current  calculation). 
The  main  CR  encircles  the  ends  of  the  bar  at  a  radius  of  about 
1 3  kpc,  with  a  pair  of  bar-driven  spiral  arms  emanating  from  the 
end  of  the  bar  outside  of  the  main  CR. 

Figure  6  uses  published  kinematic  information  to  calculate 
and  plot  several  kinematic  diagnostic  curves  for  this  galaxy, 
which  we  use  to  compare  with  the  nested  resonance  features  de¬ 
termined  using  the  phase-shift  method.  The  left  panel  shows  the 
rotation  curve  from  Regan  et  al.  (1996),  based  on  the  atomic  and 
ionized  gas  kinematics,  while  the  right  panel  shows  the  calculated 
angular  speed  —  nil,  and  the  two  pattern  speeds  determined 
by  the  phase-shift  method.  We  see  from  the  frequency  curves  that 
the  inner  ILR  (IILR)  of  the  outer  pattern  (signified  by  the  inner 


Fig.  6. — Left’.  Rotation  curve  of  NGC  1530  (Regan  et  al.  1996).  Right:  Quantities  O  {solid  curve),  Llp  {dotted  curve),  and  Q  —  n/2  {dashed  curve). 
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Fig.  7. — Comparison  of  CR  radii  {circles)  derived  from  the  phase-shift  method  with  CR  bounds  determined  by  the  TW84  method  {hatched  regions )  for  NGC  936 
{left,  Merrifield  &  Kuijken  1995)  and  NGC  4596  {right;  Gerssen  et  al.  1999).  The  phase-shift  results  are  based  on  A^-band  (2.15  gm)  images  (Laurikainen  et  al.  2005). 


intersection  of  the  £lp  of  the  outer  pattern,  i.e.,  the  lower  dotted 
line,  with  the  £2  —  k! 2  curve,  or  the  dashed  curve,  at  about  1 .3  kpc 
radius)  is  simultaneously  the  corotation  radius  of  the  inner  nested 
secondary  pattern  (signified  by  the  intersection  of  the  Qp  of  the 
inner  pattern,  i.e.,  the  upper  dotted  line,  with  the  £2  curve,  or  the 
solid  line).  This  kind  of  nested  resonance  coupling  has  been  found 
for  other  observed  galaxies  (Zhang  et  al.  1993;  Erwin  &  Sparke 
1999;  Laine  et  al.  2002)  and  in  A-body  simulations  (Schwarz  1984; 
Zhang  1998;  Rautiainen  &  Salo  1999),  and  varying  theoretical 
mechanisms  have  been  proposed  for  its  interpretation  (Contopoulos 
&  Mertzanides  1977;  Tagger  et  al.  1987;  Friedli  &  Martinet  1993). 
The  coincidence  of  nested  resonances  in  the  case  of  NGC  1530 
further  confirms  the  accuracy  of  the  phase-shift  method  in  deter¬ 
mining  the  locations  of  the  CRs. 

Also  in  Figure  6,  the  outer  intersection  of  the  lower  £lp  curve 
and  the  D  —  kJ2  curve,  which  signifies  the  outer  ILR  (OILR)  at  a 
radius  of  about  5  kpc,  coincides  with  a  depression  spot  on  the  phase- 
shift  curve  of  Figure  4  (left),  indicating  a  region  of  low  angular 
momentum  exchange  between  the  bar  density  wave  and  the  disk 
matter.  The  same  radius  also  roughly  corresponds  to  a  region  of 
low  surface  density  on  the  image  of  NGC  1 530  in  Figure  4  (right), 
where  the  central  oval  and  the  straight  bar  section  seem  to  disjoin. 

The  nested  resonances  in  the  central  regions  of  intermediate  - 
to-early-type  galaxies  are  found  to  be  important  to  the  channeling 
of  matter  to  fuel  nuclear  activity  (Zhang  et  al.  1993).  Since  these 
nested  resonances  only  form  when  the  potential  near  the  central 
region  of  a  galaxy  is  deep  enough,  the  outer  mass  accretion,  which 
helps  to  transform  a  galaxy’s  morphology  from  late  to  early 
Hubble  type  (see  §  4.3  for  further  discussions),  appears  to  serve  as 
a  preparatory  stage  for  the  fueling  of  the  active  nuclei.  The  same 
mass  accretion  process  across  a  hierarchy  of  resonance  features 
may  also  have  an  important  role  in  the  formation  of  the  tight  black 
hole  mass/bulge  mass  correlation  observed  for  many  early-type 
galaxies  (Gebhardt  et  al.  2000  and  references  therein;  Merritt  & 
Ferrarese  2001  and  references  therein). 

3.3.  Comparison  with  Corotation  Radii  Determined 
Using  the  TW84  Method 

As  noted  in  §  1 ,  the  derivation  of  galaxy  pattern  speeds  and  co¬ 
rotation  radii  has  historically  been  difficult.  The  TW84  method  is 
considered  to  be  the  most  direct  approach;  it  uses  the  continuity 


equation  and  off-nuclear  spectra  to  derive  pattern  speeds,  from 
which  a  corotation  radius  may  be  inferred  if  the  rotation  curve  is 
known.  Although  widely  applied,  the  method  has  numerous  op¬ 
erational  difficulties  that  make  it  impractical  in  some  circum¬ 
stances:  (1)  it  is  not  easily  applied  to  face-on  galaxies;  (2)  the  bar 
cannot  be  end-on  or  broadside;  (3)  the  off-nuclear  spectra  can  be 
so  time-consuming  to  obtain  that  few  galaxies  have  received  the 
telescope  time  needed;  (4)  the  continuity  equation  must  be  sat¬ 
isfied  for  the  tracers  being  used,  which  is  the  reason  the  method 
has  been  applied  mainly  to  SB0  galaxies,  leaving  spirals  largely 
unstudied;  and  (5)  its  application  is  complicated  and  model- 
dependent  if  there  is  more  than  one  pattern  present. 

In  Figure  7  we  plot  a  comparison  of  the  corotation  radii  calcu¬ 
lated  by  the  phase-shift  method  (circles)  with  the  corotation  bounds 
derived  by  the  TW84  method  (hatched  regions)  for  two  early- 
type  barred  galaxies:  the  SB(rs)0+  galaxy  NGC  936  (Merrifield 
&  Kuijken  1995)  and  the  SB(rs)0/a  galaxy  NGC  4596  (Gerssen 
et  al.  1999).  The  TW84  method  as  applied  to  these  two  galaxies 
was  used  to  determine  only  the  outer  pattern  CR  radius,  while  we 
find  inner  CR  radii  as  well.  In  each  case,  the  phase-shift  method 
gives  an  outer  CR  radius  in  good  agreement  with  the  TW84 
method. 

In  Figure  8  we  plot  a  comparison  of  corotation  radii  calculated 
by  the  phase-shift  method  (circles)  with  the  corotation  bounds 
derived  with  the  TW84  method  (hatched  regions ;  Hernandez 
et  al.  2005)  for  the  late-type  barred  galaxy  M 1 00/NGC  4321.  The 
phase-shift  results  are  calculated  using  a  Spitzer  SINGS  3.6  pm 
Legacy  image  (Kennicutt  et  al.  2003).  We  find  four  well-resolved 
corotation  radii  for  this  galaxy,  whereas  Hernandez  et  al.  (2005) 
previously  derived  three  corotation  radii  using  a  2D  velocity  field 
and  restricting  their  analyses  to  the  different  regions  assumed  for 
different  modes. 

The  CR  radius  for  the  outer  spiral  pattern  determined  by  both 
the  Hernandez  et  al.  analysis  (outer  hatched  region',  140"- 160") 
and  the  phase-shift  analysis  (outer  red  circle ;  138")  falls  in  the 
middle  of  the  outer  faint  spiral  arms,  close  to  where  the  inner, 
stronger  spiral  pattern  appears  to  truncate.  For  this  CR  the  phase- 
shift  method  and  the  TW84  method  give  results  which  seem  to 
be  in  good  agreement. 

The  innermost  corotation  determined  by  the  phase-shift  method 
(13"),  which  surrounds  the  nuclear  bar-ring  pattern,  was  not 
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Fig.  8. — Left.  Phase-shift  plot  for  Ml  00  (NGC  4321)  based  on  a  deprojected  3.6  /im  SINGS  image.  Right  Comparison  of  the  CR  radii  {circles)  derived  from  the  phase- 
shift  method,  with  the  CR  bounds  determined  by  the  TW84  method  {hatched regions)  from  Hernandez  et  al.  (2005)  for  Ml 00.  The  panel  on  the  right  is  12.8'  x  12. 8'.  These 
CR  result  comparisons  are  superimposed  on  the  3.6  \i m  Spitzer  SINGS  image  of  Ml 00  (Kennicutt  et  al.  2003). 


detected  by  the  Hernandez  et  al.  (2005)  analysis,  possibly  be¬ 
cause  of  the  lack  of  spatial  resolution  for  the  data  set  used.  For 
the  intermediate  corotations,  the  phase-shift  method  found  two 
(31"  and  59"),  which  sandwich  the  first  corotation  circle  de¬ 
termined  by  Hernandez  et  al.  (32" — 61";  innermost  hatched 
region).  It  is  at  this  point  open  to  question  whether  these  two 
closely  spaced  intermediate  crossings  indicate  two  sets  of  quasi¬ 
steady  modes  with  two  different  pattern  speeds,  or  are  caused 
by  noise  or  the  non-quasi-steady  nature  of  the  pattern.  However, 
we  found  the  middle  hatched  region  (80"- 110")  as  determined 
by  Hernandez  et  al.  to  have  no  corresponding  feature  in  the 
phase-shift  calculation.  This  middle  hatched  region  was  sug¬ 
gested  by  Hernandez  et  al.  to  be  the  corotation  of  the  bar  of 
M100. 

To  deduce  the  bar  radius,  we  used  the  SINGS  image  to  per¬ 
form  a  bar-spiral  separation  using  the  methods  described  by  Buta 
et  al.  (2005),  the  results  of  which  are  shown  with  our  four  CR  cir¬ 
cles  superimposed  in  Figure  9  ( top  panels).  From  the  separated 
bar  image,  we  estimate  a  bar  radius  ranging  from  65.8",  based  on 
the  isophote  of  maximum  ellipticity,  to  68.5",  based  on  the  faintest 
detectable  bar  isophote.  The  closest  positive-to-negative  phase 
shift  crossing  to  the  bar  radius  is  r(CR3)  =  58.1",  which  is  slightly 
inside  the  bar  ends.  We  consider  our  corotation  radius  a  more 
natural  choice  for  the  bar  than  the  Hernandez  et  al.  (2005)  deter¬ 
mination,  which  is  much  further  out  at  80"  - 1 1 0"  in  the  middle  of 
the  bright  inner  spiral,  whereas  our  CR3  circle  contains  little  of  the 
main  spiral  pattern.  Our  bar-size  determination  is  only  slightly 
longer  than  the  value  (63")  obtained  by  Knapen  et  al.  (1995). 
Given  the  intrinsic  uncertainties  in  the  bar-spiral-separation  ap¬ 
proach,  these  two  determinations  can  be  considered  in  perfect 
agreement.  Our  CR3  determination  gives  a  value  for  the  main 
bar  feature  which  is  also  close  to  the  CRbar  values  determined  by 
Arsenault  et  al.  (1988)  and  Rand  (1995). 

The  bottom  panels  of  Figure  9  present  the  original  images 
(without  bar-spiral  separation)  zoomed  in  by  a  linear  factor  of  2 
and  4,  respectively,  compared  to  the  top  panels.  From  these  two 
panels  we  can  see  more  clearly  that  the  innermost  corotation 
circle  (CR^  encloses  the  strong  secondary  bar.  Between  the  next 
two  CRs  (CR2  and  CR3)  there  appear  to  be  faint  spiral  structures. 
The  correspondence  between  the  CR  circles  and  morphological 
patterns  highlights  the  possibility  that  all  four  CRs  could  be  real, 


although  we  do  not  exclude  the  possibility  that  some  of  the 
structures  in  the  intermediate  regions  may  be  transient. 

We  mention  also  that  Oey  et  al.  (2003)  used  an  H  n  region  spa¬ 
tial  isochrone  method  to  determine  the  corotation  radii  for  Ml  00. 
For  the  outermost  corotation  they  found  a  value  (154")  similar  to 
that  of  Hernandez  et  al.  (2005)  using  the  TW84  method  (140"- 
160"),  as  well  as  ours  (138")  using  the  phase-shift  method,  but 
they  did  not  find  in  their  data  the  intermediate  corotation  location 
found  by  Hernandez  et  al.  and  others  (80"- 1 10"),  which  is  also 
absent  in  our  results.  The  H  n  region  isochrones  are  expected  to 
be  fairly  reliable  for  determining  the  corotation  crossings,  since 
the  density-wave-triggered  star  formation  behavior  has  been  well 
confirmed  from  many  previous  studies. 

In  the  process  of  calculation  for  this  galaxy  we  have  found 
that  the  large  amount  of  star  formation  in  the  Ml 00  image  does 
not  affect  the  phase-shift  results.  If  we  remove  many  of  the  star¬ 
forming  regions,  we  get  almost  the  same  CR  locations.  This  is 
apparently  related  to  the  fact  that  the  potential  as  calculated 
through  the  Poisson  integral  is  determined  by  the  global  distri¬ 
bution  of  matter  and  is  insensitive  to  the  small-scale  local  var¬ 
iations  in  intensity. 

3.4.  Leading  Arms  in  a  Grand-Design  Spiral  Galaxy 

The  identification  of  a  leading  spiral  pattern  is  equivalent  to 
determining  the  sense  of  pattern  rotation  (i.e.,  clockwise  or  coun¬ 
terclockwise,  which  is  not  the  same  as  the  sense  of  winding,  i.e., 
S-  or  Z-shaped  spirals).  Once  we  can  be  sure  which  way  the  pat¬ 
tern  is  rotating,  then  the  leading  spiral  or  bar  identification  can  be 
done  mostly  by  eye.  In  that  sense  the  leading  spiral  identification 
does  not  have  to  involve  phase-shift  calculations  at  all,  if  there  is 
some  other  means  to  identify  the  sense  of  pattern  rotation. 

In  the  course  of  our  study  we  have  discovered  an  unexpected 
application  of  the  phase-shift  approach  for  the  confirmation  of 
the  leading  nature  of  the  two  outer  spiral  arms  of  the  SA(r)a 
galaxy  NGC  4622  (Fig.  10).  From  the  sense  of  winding  of  the 
outer  arms,  it  would  initially  be  placed  in  the  “S”  category;  thus, 
we  should  choose  to  use  the  “plus”  plot  of  phase  shift  versus 
radius,  if  we  assume  that  the  galaxy  is  rotating  counterclockwise 
and  thus  the  outer  arms  are  trailing.  With  this  choice  there  are 
positive-to-negative  crossings  at  28.7"  and  47.0".  The  latter  ra¬ 
dius  occurs  in  a  no-data  area  and  so  is  ignored.  There  is  a  strong 


Fig.  9. —  Top  left  Bar-separated  SINGS  image  of  the  bright  inner  region  of  NGC  432 1  superimposed  with  the  four  corotation  circles  determined  using  the  phase-shift 
method.  The  size  of  the  image  is  6.4'  x  6.4'.  Top  right :  Spiral-separated  SINGS  image  of  NGC  4321  of  the  same  region  as  at  left  superimposed  with  the  four  corotation 
circles  determined  using  the  phase-shift  method.  The  size  of  this  image  is  6.4'  x  6.4'.  Bottom  left  SINGS  image  (without  bar-spiral  separation)  of  NGC  432 1  with  a  factor 
of  2  linear  zoom  compared  to  the  top  panels  (3 .2'  x  3 .2')  superimposed  with  the  central  three  corotation  circles  determined  using  the  phase-shift  method.  Bottom  right  SINGS 
image  (without  bar-spiral  separation)  of  NGC  4321  with  a  factor  of  4  linear  zoom  compared  to  the  top  panels  (1 .6'  x  1 .6')  superimposed  with  the  central  two  corotation  circles 
determined  using  the  phase-shift  method. 


Fig.  10. — Left  Phase  shiftvs.  radius  for  the  SA(r)a  galaxy  NGC  4622.  The  phase  shift  is  plotted  as  — </>0  because  the  galaxy  is  an  S-shaped  spiral  where  we  have  selected 
the  “minus”  plot  rather  than  the  usual  “plus”  plot.  Thus,  the  sign  convention  is  such  that  the  galactic  rotation  direction  is  clockwise.  Note  also  that  this  phase-shift  plot  was 
made  assuming  a  two-armed  spiral,  which  is  appropriate  for  the  outer  region.  For  the  inner  single-armed  pattern,  the  actual  phase-shift  value  should  be  doubled,  although 
the  zero  crossings  would  remain  the  same.  Right  Log-scale  /-band  image  of  NGC  4622  with  superimposed  corotation  circles  (red)  of  radii  24. 1"  and  35.8"  and  the  major 
negative-to-positive  crossing  (green  circle )  where  the  inner  trailing  pattern  and  the  outer  leading  pattern  decouple.  This  image  has  not  been  corrected  for  the  small 
inclination  (19°)  derived  by  Buta  et  al.  (2003). 
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negative-to-positive  crossing  in  the  outer  arms  at  35.8"  if  we 
assume  this  sense  of  galactic  rotation.  The  28.7"  circle  lies  be¬ 
tween  the  inner  and  outer  arms. 

However,  if  we  assume  that  the  inner  arm  trails,  the  outer  arms 
lead,  and  the  galaxy  rotates  clockwise,  as  was  recently  demon¬ 
strated  in  Buta  et  al.  (2003),  then  there  is  a  strong  positive-to- 
negative  crossing  at  35.8"  and  a  weaker  one  at  24. 1".  The  former 
crossing  is  almost  identical  to  the  corotation  radius  deduced  in 
Buta  et  al.  (2003)  by  examining  phase  differences  between  B- ,  V- , 
and  I- band  m  =  2  Fourier  components.  So  it  appears  that  the 
clockwise  rotation  assumption  leads  to  a  much  more  meaningful 
interpretation  of  the  phase-shift  data,  with  results  which  also 
agree  with  the  independent  analysis  based  on  the  dust /velocity 
field  information  used  by  Buta  et  al.  (2003). 

If  we  did  not  have  the  dust /velocity  field  information,  inde¬ 
pendently  we  would  not  have  been  able  to  deduce  that  the  outer 
arms  in  NGC  4622  are  leading.  Without  the  dust  lane  informa¬ 
tion  and  the  corroboration  of  the  inner  spiral  pattern,  it  would  be 
a  much  weaker  case  to  argue  that  the  outer  arms  are  leading  just 
from  the  phase-shift  results,  because  we  really  do  not  know  where 
corotation  should  be.  To  make  sense  of  the  phase-shift  plot  we 
could  also  have  chosen  the  positive-to-negative  crossing  at  28.7" 
from  the  “plus”  plot  by  assuming  counterclockwise  galactic  rotation. 

Since  we  had  to  use  the  “minus”  plot  for  an  S-shaped  spi¬ 
ral,  and  the  corotation  crossing  determined  from  a  positive-to- 
negative  transition  agrees  with  other  independent  determinations, 
that  confirms  that  we  have  identified  the  sense  of  rotation  of  the 
pattern  correctly,  since  even  for  leading  patterns,  the  angular  mo¬ 
mentum  density  of  the  pattern  is  still  negative  relative  to  the  disk 
rotation  inside  corotation  (i.e.,  the  pattern  rotates  slower  than  the 
stars  inside  corotation),  and  vice  versa  outside,  so  a  positive-to- 
negative  phase-shift  change  means  that  the  pattern  deposits  neg¬ 
ative  angular  momentum  to  the  disk  inside  corotation,  and  vice 
versa  outside,  so  the  wave  can  grow  by  an  enhanced  amplitude 
(Zhang  1998). 

This  does  not  imply  that  the  leading  spiral  of  concern  trans¬ 
fers  angular  momentum  outward.  It  implies  only  that,  since  in¬ 
side  corotation  the  wave  rotates  slower  than  the  stars  and  outside 
corotation  it  is  the  other  way  around,  if  we  take  angular  momen¬ 
tum  away  from  the  stars,  more  stars  can  join  the  wave  motion 
(and  thus  enhance  the  wave  amplitude).  For  a  positive  phase 
shift  inside  corotation,  the  wave  torques  the  stellar  density  back¬ 
ward,  so  the  wave  is  self-enhancing  for  such  a  phase-shift  dis¬ 
tribution.  It  says  nothing  about  the  radial  direction  of  angular 
momentum  flow,  only  the  angular  momentum  exchange  between 
the  wave  and  the  disk  stars  in  an  annulus.  Of  course,  the  angular 
momentum  lost  by  the  stars  at  a  particular  annulus  inside  coro¬ 
tation  is  gained  by  all  the  other  stars  which  made  up  the  wave 
motion  collectively,  since  gravity  is  long-range.  We  do  not  need 
an  actual  “flow”  to  take  away  the  angular  momentum  from  the 
stars  in  an  annulus,  although  in  the  case  of  a  true  spiral  mode 
such  a  radial  flow  of  angular  momentum  does  happen  (Zhang 
1998) — and  it  might  be  happening  here,  too,  except  that  the  flow 
direction  could  be  either  inward  or  outward.  The  important  point 
here  is  not  the  radial  angular  momentum  flow  direction  but  rather 
the  radial  gradient  of  such  a  flow  that  is  revealed  by  the  sign  of 
the  phase  shift.  If  there  is  such  a  flow,  the  gradient  of  the  flow — 
i.e.,  the  difference  between  the  entering  and  exiting  ports — must 
be  such  that  the  stars  in  an  annulus  inside  corotation  lose  angular 
momentum,  in  order  to  be  consistent  with  a  potential-density 
phase  shift  of  positive  sign. 

NGC  4622  is  a  unique  case  among  the  galaxies  we  have  an¬ 
alyzed  so  far,  in  that  the  density  wave  patterns  in  this  galaxy 
(both  the  inner  single  trailing  pattern,  the  intermediate  “arc”  just 


outside  the  CR1;  and  the  outer  two  leading  arms)  are  likely  to  be 
entirely  a  result  of  tidal  interactions.  In  a  recent  paper,  Byrd  et  al. 
(2006)  found  that  both  the  inner  and  outer  density  wave  features 
can  be  reproduced  by  the  past  plunging  passage  of  a  small  com¬ 
panion  1/100  of  the  mass  of  NGC  4622.  Furthermore,  this  model 
predicts  that  the  halo  of  this  galaxy  is  10  times  more  massive 
than  the  disk,  which  is  likely  to  be  one  of  the  reasons  that  the  disk 
intrinsic  modes  are  not  prominent.  Partly  as  a  result  of  the  mas¬ 
sive  halo,  as  well  as  the  spherical  bulge  component  in  the  central 
region,  the  phase- shift  results  for  this  galaxy  in  the  inner  region 
are  pure  noise  and  are  not  coherent  modal  responses.  Our  inner 
crossing  (24")  is  close  to  the  value  (21.5")  determined  by  other 
independent  methods  as  outlined  in  Buta  et  al.  (2003)  and  Byrd 
et  al.  (2006),  although  not  identical,  perhaps  reflecting  the  fact 
that  the  density  waves  in  this  galaxy  are  not  intrinsic  modes  of 
the  disk  and  have  not  reached  a  quasi-steady  state. 

In  general,  one  should  always  start  by  assuming  that  the  outer 
pattern  is  trailing  (which  had  been  found  to  be  the  case  for  all  spon¬ 
taneously  formed  TV-body  spiral  modes),  and  then  see  whether  the 
phase-shift  plot  and  its  predicted  corotation  radii  make  sense  when 
overlaid  on  the  image.  Of  course,  in  cases  like  NGC  4622,  where 
the  inner  and  outer  patterns  have  oppositely  winding  senses,  it  is 
obvious  that  one  set  of  patterns  must  be  leading.  The  question  then 
becomes  which  one.  This  still  assumes  that  the  inner  and  outer 
patterns  rotate  in  the  same  direction.  There  exist  cases  where 
even  the  rotation  directions  of  the  inner  and  outer  parts  of  the 
disk  stars  (and  possibly  also  density  waves)  are  different.  So,  in 
general,  we  will  need  supplementary  information,  in  addition  to 
phase-shift  calculations,  to  determine  the  sense  of  galactic  ro¬ 
tation  and  the  sense  of  leading  and  trailing  density  wave  patterns. 

An  additional  point  to  note  is  that  for  this  galaxy  the  spiral 
patterns  present  are  likely  to  be  only  transient  density  waves ,  in¬ 
stead  of  density  wave  modes.  That  the  phase-shift  method  also 
works  in  this  case  is  fortuitous  for  us,  and  in  other  cases  of  tran¬ 
sient  waves  it  should  be  expected  that  the  method  may  not  give 
correct  corotation  radii  if  the  pattern  is  going  through  rapid  mor¬ 
phological  changes. 

3.5.  Trend  of  Bar /Corotation  Sizes  versus  Morphological  Type 

An  important  issue  in  density  wave  pattern  speed  studies  is  the 
ratio  of  the  corotation  radius,  r(CR),  to  the  bar  radius,  r{ bar). 
If  this  ratio  is  less  than  1 .4,  a  bar  is  considered  “fast,”  while  if  it 
is  larger  than  1.4,  a  bar  is  considered  “slow”  (Debattista  & 
Sellwood  2000).  Here  the  words  “fast”  or  “slow”  mean  compar¬ 
ing  the  bar  pattern  speeds  to  the  local  circular  speeds  at  the  lo¬ 
cation  of  the  bar  ends. 

We  have  examined  the  ratio  of  r(CR)/r(bar)  in  the  OSUBGS 
sample  using  images  where  the  bar  has  been  separated  from  the 
spiral  using  a  systematic  approach  (Buta  et  al.  2005).  The  result 
is  given  in  Figure  1 1 .  Two  definitions  of  the  bar  radius  are  used: 
the  radius  of  maximum  ellipticity  {left)  and  the  radius  of  the 
lowest  intensity  bar  contour  {right).  The  true  bar  length  should 
be  bracketed  by  these  two  approaches.  The  plots  include  only 
those  objects  for  which  the  relative  gravitational  bar  torque  strength 
Qb  >  0.25  (Buta  et  al.  2005),  to  minimize  some  of  the  confusion 
associated  with  weaker  bars. 

From  Figure  1 1  we  see  that,  in  general,  the  bar  length  increases 
as  the  Hubble  type  of  a  galaxy  changes  from  late  to  early,  converg¬ 
ing  toward  r{ bar)  «  r(CR)  for  the  earliest  morphological  types,  a 
result  also  previously  found  by  other  investigators  (Elmegreen  & 
Elmegreen  1985;  Rautiainen  et  al.  2005;  and  references  therein). 

Among  the  galaxies  calculated  for  the  OSUBGS  sample,  we 
have  encountered  a  few  cases  in  which  the  bar  appears  to  extend 
beyond  its  CR  circle  (the  inner  bar  of  NGC  4321,  which  we  have 
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Fig.  1 1 . — Ratio  of  phase-shift-derived  CR  radius  to  bar  radius  for  36  strongly  barred  galaxies  in  the  OSUBGS,  based  on  bars  that  have  been  separated  from  their  spirals 
(Buta  et  al.  2005).  Left :  In  this  plot  the  bar  radius  was  derived  from  the  isophote  having  maximum  ellipticity.  Right:  Bar  radius  derived  from  the  faintest  distinct  isophote. 
For  this  definition,  over  the  entire  sample,  (r(CR)/r(bar))  =  1.26  ±  0.1 1.  In  these  plots,  T(  RC3)  is  the  numerically  coded  RC3  Hubble  type,  where  higher  T represents 
later  type,  i.e.,  T  =  0  is  SO/a,  T  =  1  is  Sa,  T  =  2  is  Sab,  etc. 


encountered  before,  is  a  mild  case  of  this  category),  a  phenom¬ 
enon  which  seems  to  challenge  the  so-far  well-accepted  wisdom 
of  the  extent  of  bars  in  galaxies  (Contopoulos  1980).  Figure  12 
gives  one  such  example:  the  SB(s)0+  galaxy  NGC  4665,  shown 
with  an  H- band  image  from  the  OSUBGS  sample.  From  this  fig¬ 
ure  it  is  obvious  that  this  well-defined  bar  pattern  extends  beyond 
its  own  corotation  radius. 

An  even  more  extended  bar  is  shown  in  Figure  1 3  for  OSUBGS 
galaxy  NGC  1073.  Here  the  bar  appears  to  extend  all  the  way  to 
the  negative-to-positive  crossing  on  the  phase-shift  plot,  a  location 
likely  to  coincide  with  the  OLR  of  the  inner  bar.  Other  cases 
which  appear  to  fall  into  this  category  include  OSUBGS  galaxies 
NGC  3513,  7741,  4902,  4781,  and  4579.  We  comment  more  on 
the  possible  dynamical  mechanisms  enabling  the  longer  bar  ex¬ 
tension  in  §  4.1. 

Note  that  for  the  cases  of  both  NGC  4665  and  NGC  1073, 
there  is  a  slight  ambiguity  of  how  one  defines  the  “bar.”  For 
NGC  4665,  the  bar  is  slightly  twisted  in  its  outer  region  to  appear 


somewhat  spiral-like.  For  NGC  1073,  on  the  other  hand,  at  the 
nominal  corotation  circle  there  is  a  slight  “pinch”  in  the  width  of 
the  pattern,  with  the  inner  bar  within  the  corotation  circle  having 
nearly  ellipsoidal  isophotes,  whereas  the  morphology  of  the  outer 
extension  of  the  bar  is  more  elongated  and  has  reduced  surface 
density  compared  to  the  portion  inside  r(CR).  Despite  this  am¬ 
biguity  in  classification,  one  thing  is  certain:  that  the  entire  “bar” 
structure,  both  inside  and  outside  corotation,  should  be  rotating 
with  a  single  pattern  speed. 

We  wrap  up  this  section  of  the  paper  by  summarizing  in  Table  1 
the  numerical  values  of  the  bar  lengths  and  corotation  radii  of  the 
few  galaxies  we  have  analyzed  quantitatively  in  this  work.  The 
average  of  the  ratios  of  CR  radius  to  bar  radius  for  this  small 
sample  of  seven  barred  galaxies  (excluding  NGC  4622)  are 
(r(CR)/r(bar))  =  1.01  =b  0.15  if  we  use  the  radius  of  maximum 
ellipticity  for  the  bar  and  (r(CR)/r(bar)}  =  0.87  zb  0.14  if  we 
use  the  “faintest  isophote”  for  the  bar.  The  detailed  values  for 
other  galaxies  in  the  list  of  36  strongly  barred  galaxies,  as  well  as 


Fig.  12. — Left :  Phase  shift  vs.  radius  for  the  SB(s)0+  galaxy  NGC  4665,  indicating  one  major  positive-to-negative  crossing,  as  well  as  an  additional  inner  crossing. 
Right :  Deprojected  log-scale  //-band  image  of  NGC  4665  with  superimposed  CR  circles  (red). 
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Fig.  1 3 . — Left.  Phase  shift  vs.  radius  for  the  SB(rs)cd  galaxy  NGC  1073.  Two  CR  radii  are  indicated,  neither  of  which  coincides  with  the  bar  ends  (short  horizontal  line). 
Instead,  these  ends  appear  to  be  close  to  a  major  negative-to-positive  crossing.  Right  Deprojected  log-scale  H-  band  image  of  NGC  1073  with  superimposed  CR  circles 
(red).  The  green  circle  shows  the  major  negative-to-positive  crossing. 


for  other  galaxies  in  the  OSUBGS,  will  be  given  in  the  sequel 
paper  (R.  J.  Buta  &  X.  Zhang  2007,  in  preparation). 

4.  ANALYSIS  AND  DISCUSSION 
4. 1 .  Extent  of  the  Spiral  and  Bar  Patterns 

Although  during  the  early  stage  of  the  density  wave  devel¬ 
opment  it  was  suggested  (Lin  1970)  that  the  spiral  pattern  should 
terminate  at  corotation,  this  conclusion  was  subsequently  su¬ 
perseded  by  improved  theoretical  analyses  (Mark  1976;  Lin  & 

Lau  1979;  Toomre  1981),  observational  results  (Canzian  1998  and 
references  therein),  and  TV-body  simulations  (Donner  &  Thomasson 
1994;  Zhang  1996),  which  all  showed  that  spiral  density  waves 
can  extend  to  as  far  as  the  OLR  of  the  pattern,  although  a  fraction 
of  the  observed  spiral  galaxies  indeed  displayed  evidence  of  a 
significant  reduction  of  surface  brightness  beyond  the  corotation 
radii  (Zhang  et  al.  1993;  Oey  et  al.  2003). 

TABLE  1 

Summary  of  Corotation  and  Bar  Radii 


r(CRi)  r(CR2)  r(CR3)  KCR4)  r(bar) 

Galaxy  (arcsec)  (arcsec)  (arcsec)  (arcsec)  (arcsec)  r(CRbar)/^(bar) 

(1)  (2)  (3)  (4)  (5)  (6)  (7) 


NGC  936 .  15.9  52.0*  ...  ...  47.8-55.0  1.09-0.95 

NGC  1073 .  38.0*  89.3  ...  ...  48.1-58.0  0.79-0.66 

NGC  1530 .  5.8  67.1*  ...  ...  55.1-66.1  1.22-1.02 

NGC  4314 .  9.0  76.5*  ...  ...  73.5-81.8  1.04-0.94 

NGC  4321 .  12.5  31.0  58.1*  138.0  65.8-68.5  0.88-0.85 

NGC  4596 .  13.3  69.9*  ...  ...  63.0-73.6  1.11-0.95 

NGC  4622a .  24.1  35.8 

NGC  4665 .  5.4  47.2*  ...  ...  51.0-66.0  0.92-0.71 

NGC  5247 .  13.9  86.9 


Notes. — Col.  (1):  Galaxy  name.  Cols.  (2)— (5):  Radii  of  positive-to-negative  phase-shift  crossings,  taken  to  be  CR  radii.  An  asterisk 
indicates  that  the  number  is  taken  to  be  a  bar  CR  radius.  All  of  the  galaxies  have  inner  CRs  associated  with  inner  spirals,  secondary  bars, 
or  ovals,  and  NGC  4321  shows  four  prominent  CR  crossings  corresponding  to  well-recognized  morphological  features.  Several  gal¬ 
axies  show  additional  crossings  in  the  very  centers  of  the  nuclear  region  which  were  barely  resolved  by  the  phase-shift  calculations 
(NGC  4314,  NGC  5247,  NGC  1530,  NGC  4321,  and  NGC  1073);  the  values  for  these  crossings  were  not  explicitly  given  due  to  the 
insufficient  resolution  at  these  galactic  central  regions.  Col.  (6):  Estimates  of  the  bar  radius  based  on  two  methods:  fitted  radius  of  the 
ellipse  of  maximum  ellipticity  (smaller  value)  and  radius  of  the  faintest  detectable  isophote  (larger  value),  based  on  an  image  where 
the  bar  has  been  separated  from  its  associated  spirals  or  rings.  Several  galaxies  have  a  certain  degree  of  asymmetry  between  the  two 
sides  of  the  bar,  such  as  NGC  1073.  Col.  (7):  Ratio  of  CR  to  bar  radius  based  on  the  range  in  col.  (6). 
a  The  CR  radii  are  based  on  a  graph  of  —  (f)0  vs.  radius  (see  text). 


For  barred  galaxies,  a  similar  conclusion  (that  the  pattern  ter¬ 
minates  at  its  own  corotation)  was  reached  by  Contopoulos  (1980), 
based  on  orbit  analysis  under  an  applied  potential  to  construct  “self- 
consistent”  models  of  barred  galaxies.  It  was  found  that  the 
phases  of  the  periodic  orbits  thus  obtained  usually  changed  rap¬ 
idly  when  approaching  the  corotation  radius  (the  phase  change 
starts  to  be  significant  around  the  4 : 1  resonance,  which  lies  within 
corotation),  making  such  self-consistent  model  constructions  dif¬ 
ficult.  This  was  the  underlying  reason  for  Contopoulos’s  (1980) 
conclusion  that  a  bar  should  terminate  at  corotation. 

Since  it  is  by  now  well  established  that  a  spiral  pattern  can 
extend  beyond  its  corotation  radius,  and  since  bars  are  density 
wave  patterns  just  like  spirals,  apart  from  the  fact  that  they  reside 
in  the  central  regions  of  galaxies,  one  naturally  suspects  that  a 
bar  might  be  able  to  extend  beyond  its  own  corotation  radius  as 
well.  In  fact,  this  possibility  was  hinted  at  already  in  Contopoulos 
(1980),  where  the  orbital  behaviors  of  stars  near  the  ILRs  of 
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galaxies  were  analyzed  in  addition  to  those  near  the  corotation. 
In  particular,  for  galaxies  which  possess  two  ILRs,  even  though 
the  inner  bar  normally  ends  within  the  IILR,  in  certain  cases  “the 
inner  bar  can  extend  up  to  the  outer  ILR  if  its  amplitude  drops 
sufficiently  fast”  (Contopoulos  1980). 

For  nested  resonances,  it  often  happens  that  the  corotation 
of  the  inner  secondary  pattern  is  close  to  the  location  of  the  IILR 
of  the  outer  pattern  (as  we  have  seen  for  the  case  of  NGC  1530) 
and  the  OLR  of  the  secondary  pattern  is  close  to  the  OILR  of  the 
outer  pattern  (for  the  case  of  NGC  1530,  although  the  OLR  of 
the  inner  pattern  resides  within  the  OILR  of  the  outer  pattern). 
Therefore,  Contopoulos’s  conclusion  that  the  secondary  pattern 
can  sometimes  fill  in  the  region  between  two  ILRs  implicitly  ad¬ 
mits  that  the  inner  secondary  bar  can  extend  beyond  its  own  co¬ 
rotation  and  reach  all  the  way  to  its  own  OLR,  and  perhaps  even 
beyond  (if  the  OLR  of  the  inner  pattern  lies  within  the  OILR  of 
the  outer  pattern). 

The  evidence  that  bars  can  extend  beyond  their  corotation  has 
in  fact  already  shown  up  in  earlier  self-consistent  A-body  sim¬ 
ulations.  Sparke  &  Sellwood  (1987)  produced  one  of  the  earliest 
self-sustained  2D  A-body  bar  modes,  which  settled  into  a  quasi¬ 
steady  state  after  first  emerging  as  a  spiral  pattern.  During  much 
of  the  simulation  run  after  the  mode  reached  a  nonlinear  satu¬ 
ration  amplitude,  the  extent  of  the  bar  was,  in  fact,  longer  than 
the  corotation  radius  at  the  corresponding  time  step.  For  exam¬ 
ple,  from  their  Figure  3  a ,  we  can  see  that  the  inner  bar  mode  has 
contours  extending  to  about  1 .7  in  radius  at  time  step  40,  whereas 
the  corotation  radius  for  the  same  time  step  is  1 .  In  a  later  stage 
(step  103)  of  the  simulation  (their  Fig.  3b)  the  corotation  grows 
to  1 .3  (due  to  the  slowing  down  of  the  bar  pattern);  the  bar,  on  the 
other  hand,  has  reduced  in  length,  so  it  is  indeed  enclosed  within 
the  corotation  radius  at  this  late  stage  of  the  simulation.  Note  that 
even  though  the  morphology  of  the  bar  and  spiral  pattern  in  the 
Sparke  &  Sellwood  (1987)  simulation  is  evolving  continuously, 
this  is  the  evolution  of  the  quasi-steady  modes  as  a  result  of 
the  evolution  of  the  underlying  basic-state  distribution  due  to  the 
wave-induced  angular  momentum  transport  and  exchange.  The 
wave  modes  themselves  achieved  a  quasi-steady  state  through¬ 
out  most  of  the  simulation  run,  as  reflected  in  the  well-defined 
pattern  speeds  that  precisely  delineate  the  extents  of  the  inner 
and  outer  modes  (see  their  Fig.  3).  That  the  Sparke  &  Sellwood 
(1987)  simulation  had  the  longest  bar  length  (compared  to  the 
size  of  the  modal  corotation  radius)  at  the  intermediate  stage  of 
the  simulation  is  also  consistent  with  our  finding  that  the  small¬ 
est  r(CR)/r(bar)  ratios  are  observed  for  intermediate  Hubble 
type  galaxies. 

We  further  argue  in  §  4.2  that  the  fact  that  bars  and  spirals  can 
extend  beyond  their  own  corotation  radii  is  partly  due  to  the  fact 
that  when  self-consistent  patterns  and  collective  effects  are  taken 
into  account,  many  of  the  conclusions  reached  through  passive  or¬ 
bit  calculations  are  no  longer  valid.  As  a  matter  of  fact,  in  general 
there  will  not  be  a  clear  dividing  line  between  bars  and  spirals 
(apart  from  the  fact  that  bars  tend  to  occur  in  the  inner  regions  of 
galaxies  and  spirals  in  the  outer  regions);  most  realistic  bars  in 
physical  galaxies  (including  the  more  than  1 00  galaxies  we  have 
examined  so  far  in  the  OSUBGS  sample)  all  possess  a  certain 
degree  of  skewness,3  which  enables  the  presence  of  the  phase 
shift,  and  the  associated  collective  dissipation  effect,  which  leads 


3  Certain  very  early  type  galaxies  may  contain  a  very  limited  amount  of 
skewness  (and  thus  a  very  small  resulting  phase  shift),  such  as  many  SO  and  disky 
elliptical  galaxies.  These  galaxies  can  be  considered  old  or  even  close  to  being 
“dead”  in  the  evolutionary  sense,  in  that  the  phase-shift-induced  secular  evo¬ 
lution  has  practically  halted. 


to  the  secular  evolution  of  the  basic  state  (see  §  4.3).  In  cases  of 
bar-driven  spirals  it  might  even  be  hard  to  tell  where  the  bar  ends 
and  the  spiral  starts  (such  as  in  NGC  4321),  especially  during  the 
intermediate  stage  of  the  evolution  when  the  spiral  is  being  trans¬ 
formed  into  a  bar-driven  spiral.  For  other  galaxies  (such  as  NGC 
1073,  NGC  4665,  the  inner  mode  of  NGC  1530,  and  the  simu¬ 
lated  galaxy  disk  in  Sparke  &  Sellwood  [1987]),  the  scenario  of 
the  bar  extending  beyond  corotation  is  achieved  through  the  de¬ 
coupling  of  the  inner  and  outer  nested  modes. 

4.2.  Orbits  as  Building  Blocks  ”  of  Large-Scale 
Density  Wave  Patterns 

Stellar  trajectories  in  a  self-consistent  spiral  pattern  can  show 
drastically  different  behavior  than  periodic  or  quasi-periodic 
orbits  calculated  in  an  applied  spiral  potential.  While  both  groups 
become  chaotic  at  large  spiral  amplitudes,  the  former  also  exhibit 
secular  mean-radius  change,  which  is  not  displayed  in  the  latter. 
Such  a  difference  in  behavior  is  due  to  the  constraint  of  the 
potential-density  phase  shift  applied  by  the  self-consistent  spiral 
pattern  on  the  individual  stars’  trajectories,  which  causes  the  spi¬ 
ral  wave  to  steepen  into  large-scale  (collisionless)  spiral  shocks. 
This  further  leads  to  the  irreversible  conversion  of  stellar  orbital 
motion  energy  into  random  motion  energy,  as  well  as  into  the  en¬ 
ergy  of  the  wave,  through  a  local  gravitational  instability  of  the 
streaming  disk  material  at  the  spiral  wave  crest  (Zhang  1 996, 1 998). 
We  now  present  a  quantitative  example  to  illustrate  the  difference 
in  behavior  between  a  passive  orbit  calculated  in  a  rigidly  rotating, 
applied  spiral /bar  potential  and  a  stellar  trajectory  obtained  from  a 
self-consistent  A-body  calculation. 

In  Figure  14  we  plot  the  nonlinear  orbit  solution,  together 
with  the  corresponding  surface  of  section,  for  spiral  forcing  of 
strength  10%,  20%,  and  30%,  under  an  axisymmetric  mean  po¬ 
tential  similar  to  that  in  the  solar  neighborhood.  The  equations 
used  for  calculating  the  orbits  are  the  standard  nonlinear  coupled 
second-order  ordinary  differential  equations  in  the  corotating 
frame  of  the  spiral  (Binney  &  Tremaine  1987,  eqs.  [3- 107a]  and 
[93-1 07b]),  i.e., 

f  -  rtj)2  =  -  +  2 r<j>Qp  +  f l2pr,  (6) 

1 

ref)  +  2r<fi  = - —  —  2rttp,  (7) 

r  ocp 

and  with  the  forcing  potential  being  of  the  form 

(2  In 

2</>-f  — — 7j,  (8) 

where  vc  is  the  circular  speed,  Qp  is  the  pattern  speed,  and  i  is 
the  pitch  angle  of  the  spiral.  In  all  the  calculations  performed  for 
Figure  14,  we  have  used  r  =  8.5  kpc,  vc  =  220  km  s-1,  Dp  = 
13.5  km  s-1  kpc-1,  and  i  =  20°.  The  fractional  forcing  strength / 
is  defined  as /  =  2v4/[tan  i(v%)\. 

We  observe  from  Figure  14  that  for  all  the  orbits  calculated, 
there  is  no  secular  change  in  the  mean  orbital  radius,  which  is  a 
result  consistent  with  the  constancy  of  the  Jacobi  integral  in  a  rigidly 
rotating  spiral  potential  (Binney  &  Tremaine  1987,  eq.  [3-88]).  Fur¬ 
thermore,  we  notice  that  orbits  are  only  regular  at  small  spiral¬ 
forcing  amplitude.  For  forcing  amplitudes  which  exceed  20%,  the 
orbits  gradually  become  chaotic,  and  their  traces  on  the  surface 
of  section  begin  to  cover  an  entire  2D  area. 
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Fig.  14. — Orbit  solution,  plotted  in  the  pattern-stationary  frame,  and  surface  of  section  for  spiral  forcing  strengths  of  10%  (a,  b ),  20%  (c,  d ),  and  30%  ( e,f ).  Parameters 


common  to  these  three  cases  are  the  pitch  angle  of  the  spiral  i  —  20° ,  Qp  —  13.5  km  s 
The  test  particle  was  initially  on  a  circular  orbit. 

It  is  obvious  that  a  coherent  spiral  density  distribution  cannot 
be  obtained  by  the  straightforward  superposition  of  orbits  of  the 
type  given  by  Figure  I4e.  However,  20%-30%  is  a  realistic  forc¬ 
ing  level  in  the  observed  spiral  and  barred  galaxies.  Therefore,  at 
least  in  the  case  of  (but  not  limited  to)  large-amplitude  spirals,  or- 


1  kpc  f  r(t  =  0)  =  8.5  kpc,  vc  —  220  km  s  ',  and  total  integration  time  t  —  80  Gyr. 


bits  can  no  longer  be  considered  the  “building  blocks”  of  the  global 
spiral  pattern.  Certain  collective  effects  have  to  be  present  in  order 
to  “collapse  the  chaos”  inherent  in  the  stellar  trajectories  and  arrive 
at  a  coherent  spiral  density  response.  In  later  years  Contopoulos  and 
collaborators  also  included  more  of  the  role  of  chaotic  orbits  in 
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Fig.  15. — Time  evolution  of  the  trajectory  of  a  typical  star  inside  corotation,  for  a  spontaneously  formed  7V-body  spiral  mode  simulated  in  Zhang  (1996).  The  trajectory 
sections  are  plotted  with  respect  to  the  pattern-stationary  (corotating)  frame.  The  forcing  amplitude  is  around  15%. 


constructing  galaxy  models  (Contopoulos  &  Voglis  1996  and 
references  therein;  Voglis  et  al.  2006),  although  the  full  collec¬ 
tive  effects  are  still  missing  in  these  passive  calculations.4 


4  The  chaos  present  in  the  passive  orbit  response  is  generally  of  a  different  nature 
than  the  chaotic  behavior  observed  in  collective  systems;  i.e.,  it  lacks  the  hidden  col¬ 
lective  correlations  among  the  TV-particles’  trajectories  in  collective  systems,  yet  it  dis¬ 
plays  more  coherent  patterns  and  correlation  with  respect  to  the  mean  forcing  potential. 


The  A-body  simulations  of  spontaneously  formed  spiral  modes 
reveal  that  typical  particle  trajectories  in  such  finite-amplitude, 
self-sustained  spiral  patterns  are  again  chaotic  in  the  corotating 
frame  of  the  pattern;  however,  a  new  phenomenon  not  present  in 
the  passive  orbit  calculations  is  now  also  evident.  A  typical  tra¬ 
jectory  of  a  particle  in  the  simulation  of  Zhang  (1996)  in  the  co¬ 
rotating  frame  of  a  quasi-stationary  spiral  pattern,  for  a  radial 
location  inside  corotation,  is  plotted  in  Figure  15.  This  simulation 
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run  obtained  a  spontaneously  formed  spiral  mode  of  equilibrium 
amplitude  around  15%.  It  is  clear  that  this  particle  trajectory 
covers  a  2D  area  in  the  reference  frame  corotating  with  the  den¬ 
sity  wave  pattern,  signaling  the  presence  of  chaos  (for  a  15%  forc¬ 
ing  without  chaos,  the  orbit  should  look  like  the  tubes  in  Figs.  14a 
and  14c).  The  invasion  of  chaos  is  thus  found  to  occur  at  a  much 
smaller  spiral-forcing  amplitude  for  the  spontaneously  formed 
A-body  spirals  than  for  the  non-self-consistent  case,  a  result  con¬ 
sistent  with  the  theoretical  finding  that  collective  relaxation  induced 
by  large-scale  global  instability  patterns  generally  dominates  the 
relaxation  behavior  of  stellar  trajectories  in  galaxies  (Weinberg 
1993;  Pfenniger  1986;  Zhang  1996). 

Another  new  phenomenon  evident  from  Figure  1 5  is  that  there 
is  now  a  secular  decay  of  the  mean  orbital  radius  for  stars  inside 
corotation.  A  corresponding  secular  increase  of  the  mean  orbital 
radius  for  stars  outside  corotation  is  also  found  (Zhang  1996). 
Both  of  these  secular  changes  in  mean  trajectory  sizes  are  absent 
from  the  non-self-consistent  forced  orbit  calculations. 

The  presence  of  stochasticity  in  an  individual  star’s  orbit  is  an 
integral  element  in  a  self-organized  and  self-sustained  spiral  pat¬ 
tern.  The  stochasticity  derives  from  the  exponential  sensitivity  of 
the  stellar  trajectory  to  the  small  perturbations  in  initial /boundary 
conditions,  with  the  sensitivity  itself  a  result  of  the  presence  of 
global  instabilities,  which  create  an  effective  “collisional  relaxa¬ 
tion”  condition  (Zhang  1996).  The  large-scale  coherence  of  a  spon¬ 
taneously  formed  spiral  pattern  is  nonetheless  not  sensitive  to  the 
presence  of  the  chaotic  behavior  in  an  individual  star’s  trajec¬ 
tory.  This  is  because  the  spontaneous  emergence,  and  the  sub¬ 
sequent  maintenance  and  secular  evolution,  of  these  coherent 
patterns  is  a  natural  sequence  which  corresponds  to  a  large  phase- 
space  volume  in  all  the  possible  initial  conditions  that  a  protoga- 
lactic  disk  could  choose  from.  The  large  increase  in  coarse-grained 
entropy  when  a  disk  develops  a  spiral  mode  means  that  the  pattern 
itself  is  an  attractor  to  the  evolution  of  these  particles’  collective 
motions. 

Therefore,  in  spiral  disks  we  witness  three  levels  of  organiza¬ 
tion:  deterministic  in  the  governing  equations,  chaos  in  the  indi¬ 
vidual  star’s  trajectory,  and  deterministic  again  in  the  large-scale 
patterns  formed  in  such  a  many-body  system.  Collective  dissipa¬ 
tion  is  what  makes  possible  the  existence  of  large-scale  coher¬ 
ence  in  a  sea  of  underlying  chaotic  orbits.  At  each  crossing  of  the 
spiral  or  bar  crest,  the  orbital  parameters  of  the  individual  stellar 
trajectories  are  reorganized  to  support  the  wave  pattern,  and  they 
are  completely  decorrelated  from  their  parameters  in  the  previous 
cycle  of  arm  crossing,  due  to  the  presence  of  gravitational  in¬ 
stability  at  the  arms.  The  presence  of  global  resonances  (such  as 
ILRs  and  CRs)  in  such  self-consistent  patterns  will  no  longer  act 
on  individual  orbits  but  rather  on  the  collective  motions  of  all  the 
particles  constituting  the  density  wave.5  In  a  sense,  the  collection 
of  stars  behaves  rather  like  molecules  in  a  fluid,  and  the  collisional 
shock  in  the  spiral  arms  and  bars  creates  conditions  of  particle 
interaction  akin  to  true  collisions  in  a  true  fluid,  albeit  here  not  in 
the  sense  of  two-body  collisions  but  rather  the  scatterings  of  par¬ 
ticles  in  the  collective  instability  at  the  density  wave  crest. 

So  far,  all  the  existing  orbital  resonance  calculations  and  gal¬ 
axy  model  constructions  have  been  done  under  an  applied  potential, 


5  The  fact  that  an  individual  star’s  resonance  condition  is  totally  destroyed  in  a 
disk  galaxy  containing  a  global  instability  can  be  understood  as  due  to  the  decor¬ 
relation  of  the  phases  of  the  trajectory  each  time  it  experiences  a  gravitational- 
instability-induced  scattering/collision  event  at  the  spiral  arm  crest.  Like  the 
example  of  a  child’s  swing,  both  the  exact  frequency  condition  and  the  exact  phase 
condition  need  to  be  satisfied  for  the  swing  to  gain  in  amplitude.  In  the  case  of 
individual  trajectories  in  galaxies,  it  is  the  phase  coherence  condition  that  is  being 
destroyed  by  the  density-wave-induced  gravitational  instability. 


not  the  kind  of  self-consistent  potential  in  many-particle  sys¬ 
tems.  Even  for  several  supposedly  self-consistent  calculations  us¬ 
ing  A-body  simulations,  a  closer  examination  of  the  ways  some 
of  the  final  analyses  were  conducted  revealed  that  it  was  only  the 
forced/passive  orbital  behavior  that  had  been  uncovered  in  these 
analyses.  For  example,  Athanassoula  (2003)  studied  the  reso¬ 
nance  exchanges  between  bars  and  active  halos.  Even  though  she 
had  started  out  with  a  self-consistent  three-dimensional  A-body 
simulation,  when  the  time  came  to  calculate  the  resonant  ex¬ 
changes  between  the  bar  stars  and  halo  stars,  she  had  to  halt  the 
A-body  simulation  and  revert  to  using  a  rigid ,  uniformly  rotating , 
and  enforced  bar  potential  which  best  approximated  the  bar  po¬ 
tential  she  had  obtained  at  the  end  of  the  A-body  simulation,  so  as 
to  obtain  the  sought-after  orbital  resonance  exchange  with  the  halo. 
It  is  not  surprising  that  such  a  practice  should  be  needed;  resonance 
conditions  exist  for  individual  orbits  only  under  an  applied  and 
rigidly  rotating  potential,  not  for  a  self-consistent  and  (micro¬ 
scopically)  time-fluctuating  potential  from  the  actual  A-body  sim¬ 
ulations  where  collective  effects  completely  wipe  out  the  individual 
orbit’s  resonance  behavior.  Likewise,  other  investigators  who  con¬ 
ducted  periodic  and  quasi-periodic  orbit  searches  in  self-consistent 
A-body  models  also  used  either  the  A-body  potential  at  one  instant 
of  the  simulation  (Pfenniger  &  Friedli  1991)  or  a  time-averaged 
potential  (Sparke  &  Sellwood  1987)  to  search  for  periodic  or¬ 
bit  families,  which  effectively  are  equivalent  to  using  a  forced, 
rigidly  rotating,  and  non-self-consistent  potential  over  the  orbital 
period  under  study.  It  is  no  wonder  that  the  secular  mean  trajec¬ 
tory  changes  we  have  found  in  the  true,  fully  self-consistent 
A-body  simulations  were  never  found  in  those  orbit  searches  which 
used  the  A-body  potential  at  a  fixed  instant  to  approximate  the  po¬ 
tential  over  the  entire  orbital  period. 

The  true  collective  behavior  is  encoded  in  and  enabled  by  the 
spatial  and  temporal  correlations  of  the  A-particles’  positions 
and  velocities.  While  the  phase  correlation  of  a  single  star’s  tra¬ 
jectory  with  the  mean  forcing  potential  has  been  destroyed  by  the 
presence  of  collective  instabilities,6  the  multifaceted  correlations 
among  the  entire  collection  of  A-particles  are  being  established 
to  drive  the  emergence,  maintenance,  and  long-term  evolution  of 
the  collective  instability.  An  analogy  from  a  different  field  is  that 
of  the  growth  of  a  living  organism.  Most  of  us  would  agree  that 
if  we  simply  assemble  together  the  material  content  of  a  living 
organism,  no  matter  how  closely  we  can  approximate  the  actual 
material  makeup  of  the  organism,  we  will  never  be  able  to  “breathe 
life”  into  such  an  assembly  of  parts.  The  “life  force”  of  a  living 
being  is  derived  from  the  very  self-organization  process  (i.e.,  the 
growth  process  of  the  biological  form),  through  the  multitude  of 
correlations  established  during  the  “co-evo”  process.  Similarly, 
for  galaxies  possessing  a  large-scale  density  wave  mode,  the  ac¬ 
tual  particles’  trajectories  supporting  the  pattern  acquire  global 
coherence  during  the  spontaneous  growth  process  of  the  mode, 
despite  their  chaotic  appearance  when  each  is  viewed  individu¬ 
ally.  In  some  sense,  the  apparently  “ordered”  forms  of  the  peri¬ 
odic  orbits  for  galactic  systems  are  static  and  inert,  and  as  a  result 
they  lack  the  very  life  force  (i.e.,  the  ability  to  maintain  themselves 
and  simultaneously  evolve  through  the  enabling  roles  of  self¬ 
organization  and  collective  dissipation)  we  are  talking  about, 
whereas  the  apparently  more  “chaotic-looking”  trajectories  in 


6  Note  that  this  appears  to  have  violated  the  conclusion  of  the  celebrated 
KAM  theorem  (after  Kolmogorov,  Arnold,  &  Moser),  which  prescribes  the  con¬ 
ditions  for  the  stability  of  periodic  and  quasi-periodic  orbits  under  small  pertur¬ 
bations.  We  point  out  that  the  KAM  theorem  was  derived  for  orbit  structures 
under  an  applied  potential,  not  for  self-organized  collective  instabilities.  Or,  in 
other  words,  for  self-organized  collective  processes,  the  perturbations  cannot  be 
regarded  as  small  from  the  outset  of  the  growth  of  the  instability. 
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the  A-body  systems  in  fact  contain  innate  order  and  capacity  for 
slow  transformation  derived  from  the  self-organization  process, 
which  facilitates  the  further  maintenance  and  evolution  of  the 
parent  galactic  system.  The  equilibrium  in  physical  galaxies  is  a 
dynamical  equilibrium,  maintained  as  a  balance  of  the  growth 
and  dissipation  tendencies  of  the  wave  mode  (Zhang  1 998),  rather 
than  the  static  equilibrium  described  by  the  passive  orbit  analysis 
and  model-building  approaches. 

The  analogies  we  make  above  between  a  biological  system  pos¬ 
sessing  living  organisms  and  a  galactic  system  possessing  den¬ 
sity  wave  modes  are  in  fact  not  accidental;  both  are  instances  of 
“dissipative  structures”  governed  by  the  laws  of  self-organization 
for  systems  that  are  open  and  far  from  equilibrium  (Zhang  1998; 
Glansdorff  &  Prigogine  1971). 

4.3.  Implications  for  the  Kinematic  and  Dynamical 
States  of  Nearby  Galaxies 

Spiral  and  barred  galaxies  are  expected  to  be  systems  in  near 
dynamical  equilibrium,  and  thus  their  kinematic  and  dynami¬ 
cal  states  are  partly  reflected  in  their  morphological  appearances; 
this  connection  between  morphology  and  kinematics  for  quasi- 
equilibrium  galaxies  is  the  foundation  of  our  approach  for  coro¬ 
tation  determination. 

So  far,  we  have  detected  the  phase  shifts  and  derived  corota¬ 
tion  radii  using  H- band  (1 .65  pm)  images  for  more  than  100  gal¬ 
axies  in  the  OSUBGS  sample,  which  includes  spiral  and  barred 
galaxies  of  all  morphological  types,  as  well  as  using  several  gal¬ 
axy  images  from  other  ground-based  and  space-based  data  sets. 
Our  initial  success  in  applying  this  method,  as  reflected  in  the 
observed  correspondence  between  the  corotation  radii  determined 
by  the  phase-shift  method  and  the  resonance  features  present  in 
the  original  images,  and  the  good  agreement  between  our  results 
and  many  of  the  previously  published  ones  using  other  methods, 
indirectly  confirms  that  a  large  fraction  of  the  observed  disk  gal¬ 
axies  in  the  nearby  universe  have  acquired  a  quasi- steady  state 
with  respect  to  the  formation  of  spiral/bar  modes,  since  only  in  the 
quasi- steady  state  of  the  wave  mode  can  one  expect  a  correlation 
between  the  kinematic  characteristics  (such  as  corotation  radius) 
and  the  pure  morphological  appearance  (i.e.,  the  NIR  image). 

Over  the  four  decades  since  the  advent  of  density  wave  theory, 
the  issue  of  whether  the  observed  density  wave  patterns  in  gal¬ 
axies  are  quasi-stationary  or  transient  has  always  been  a  matter 
of  debate  (see,  e.g.,  Binney  &  Tremaine  1987,  §  6.4  and  references 
therein;  Lucentini  2002).  Among  the  A-body-simulated  galaxies, 
results  also  differ  between  the  quasi-steady  patterns  (Donner  & 
Thomasson  1994;  Zhang  1998)  and  the  more  transient  patterns 
obtained  (Sellwood  &  Carlberg  1984;  Sellwood  &  Binney  2002). 
We  point  out  here  that  at  least  one  of  the  known  causes  of  the 
differing  results  in  these  A-body  simulations  is  the  differing  basic 
states  of  the  galactic  disks  being  used.  When  simulating  spiral 
patterns,  Sellwood  &  Evans  (2001)  and  Sellwood  &  Binney  (2002) 
used  a  disk  model  that  is  overstable  with  respect  to  the  formation 
of  large-scale  spiral  modes.  Therefore,  any  spiral  perturbations 
will  tend  to  die  out,  which  is  the  reason  they  only  obtained  transient 
spiral  patterns  in  these  studies.  However,  growing  observational 
evidence,  including  our  current  work,  shows  that  the  large-scale 
spiral  and  barred  patterns  in  physical  galaxies  are  in  fact  grav¬ 
itationally  unstable  global  modes ,  not  the  transient  structures  stud¬ 
ied  in  the  above  work,  which  used  overstable  basic  states.  As  a 
matter  of  fact,  even  in  some  of  the  work  conducted  by  these 
authors  themselves  (i.e.,  Sparke  &  Sellwood  1987;  Debattista  & 
Sellwood  2000),  when  a  basic-state  model  that  allowed  a  genuine 
global  mode  to  emerge  was  used  (in  this  case  mostly  bar  modes 
and  bar-driven  spirals),  they  had  invariably  observed  the  same 
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Fig.  1 6. — Radial  mass  accretion/excretion  rate  calculated  for  NGC  1530  from 
a  A -band  image  (§  3.2). 


kind  of  collective  effects  (including  secular  orbital  change)  and 
long-lasting  patterns. 

One  other  consequence  of  the  presence  of  a  phase  shift,  as  we 
have  already  seen,  is  the  secular  evolution  of  the  basic-state  mass 
distribution  for  a  galaxy  containing  large-scale,  skewed  density 
wave  patterns.  The  secular  morphological  transformation  mech¬ 
anisms  proposed  in  the  past  involve  mostly  gas  accretion  under  a 
barred  potential  (Combes  &  Sanders  1981;  Kormendy  &  Kennicutt 
2004  and  references  therein),  which  were  found  to  be  inadequate  to 
transform  the  Hubble  types  of  galaxies,  except  for  the  very  latest 
types  (Sc,  Sd),  due  to  the  paucity  of  gas  in  most  disk  galaxies. 
When  the  role  of  stars  is  also  considered  (Zhang  1996,  1998, 
1999),  a  more  coherent  and  theoretically  sound  basis  for  the  sec¬ 
ular  evolution  paradigm  emerges  that  is  powerful  enough  to  ac¬ 
count  for  both  the  morphological  transformation  of  all  galaxy 
types  along  the  Hubble  sequence  and  the  differing  evolution  rate 
observed  for  galaxies  residing  in  the  field  and  the  cluster  environ¬ 
ments  (cluster  galaxies  evolve  faster  because  of  the  large-amplitude, 
open  spiral  patterns  excited  during  the  tidal  interactions  with  neigh¬ 
boring  galaxies  and  with  the  cluster  potential,  and  the  secular  evo¬ 
lution  rate  is  effectively  proportional  to  the  wave  amplitude 
squared  and  pattern  pitch  angle  squared  [Zhang  1998,  1999]). 

In  general,  the  radial  mass  accretion /excretion  rate  across  any 
galactic  radius  can  be  calculated  from  (Zhang  1998) 


dM  x  R  f 

— (r)  = - / 

dt  vc  Jo 


2"  ,  ,.dV(r,£) 

<j>) — ^ —  d<t>, 


(9) 


where  the  surface  density  X  and  potential  V  have  been  written  in 
their  nonperturbative  form  (no  subscript  1,  as  in  the  previous 
equations),  since  in  this  calculation  we  can  use  the  overdensity 
and  potential  obtained  directly  from  NIR  images.  But  when  the 
actual  differentiation  and  integration  are  carried  out,  only  the  per¬ 
turbative  quantities  contribute  to  the  final  value  obtained  (Zhang 
1996). 

In  Figure  16  we  give  an  example  of  such  a  calculation  for  the 
SBb  galaxy  NGC  1530,  using  the  surface  density  map  for  NGC 
1530  described  in  §  3.2.  The  methodologies  needed  for  the  cal¬ 
culation  of  phase  shift  and  accretion  rates  are  different.  The  phase 
shift,  being  dimensionless,  can  be  determined  without  knowledge 
of  the  distance  and  with  reasonable  accuracy  under  the  assump¬ 
tion  of  a  constant  MIL.  However,  mass  accretion  rates  require  the 
knowledge  of  real  M/Zs,  galaxy  rotation  velocities,  and  absolute 
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lengths.  Note  that  the  shape  of  the  phase-shift  distribution  curve 
(Fig.  4,  left)  is  different  from  the  mass  accretion  rate  distribution 
(Fig.  16).  This  is  because  the  mass  accretion  rate  depends  not  only 
on  the  phase  shift  but  on  the  surface  mass  density  distribution  and 
the  distribution  of  the  density  wave  perturbation  amplitude,  as 
well  as  on  the  galaxy  rotation  curve  distribution.  The  peak  ac¬ 
cretion  rate  for  this  galaxy  is  close  to  140  M0  yr-1  near  the  end  of 
the  main  bar.  The  inner  accretion  rate  we  obtained  for  this  galaxy 
(5-10  M0  yr-1)  is  consistent  with  the  gas  accretion  rate  esti¬ 
mated  by  Regan  et  al.  ( 1 997 ;  see  their  F ig.  2),  which  is  in  the  range 
of  2-5  M0  yr-1  for  the  inner  4  kpc,  since  our  result  includes  both 
stellar  and  gaseous  accretions.  Also  note  that  our  results  indicate 
both  accretion  and  excretion  in  the  central  region  due  to  nested 
resonance  structures,  with  the  exact  locations  of  these  resonances 
expected  to  evolve  over  a  long  time  as  the  basic  state  itself  evolves. 
The  large  excretion  in  the  outer  disk  is  expected  to  contribute  to  the 
buildup  of  an  extended  outer  envelope  usually  observed  for  disk 
galaxies  and  disky  ellipticals. 

The  result  of  the  mass  accretion  calculation  above  shows  that 
for  at  least  a  fraction  of  galaxies  such  as  NGC  1530  which  have 
extremely  strong  and  open  density  wave  patterns,  the  secular  evo¬ 
lution  can  be  very  rapid  (i.e.,  a  mass  accretion  rate  of  1010  M0  per 
108  yr  is  possible),  compared  to  a  typical  galaxy  like  our  own 
Milky  Way,  where  it  takes  about  1 010  yr  to  accumulate  6  x  1 09  M0 
through  the  same  mechanism  (Zhang  1999),  since  the  Galactic 
spiral  is  both  less  intense  and  more  tightly  wound.  We  have  so  far 
calculated  the  mass  accretion  rates  for  several  other  galaxies  as 
well  and  found  that  a  peak  accretion  rate  of  several  to  several  tens 
of  solar  masses  per  year  is  typical  of  galaxies  of  intermediate 
Hubble  types.  An  associated  heating  effect  as  a  by-product  of  the 
energy  and  angular  momentum  exchange  process  between  the 
density  wave  and  the  basic-state  stars  also  results  in  the  secular 
heating  of  the  disk  stars  in  addition  to  their  secular  orbit  migration 
(Zhang  1999),  and  this  results  in  the  vertical  migration  of  the 
stellar  trajectories,  as  well  as  the  radial  migration,  which  together 
lead  to  the  building  up  of  the  bulge  star  populations  (Zhang  2003). 
We  note  here  that  some  of  the  past-proposed  bulge-building  mech¬ 
anisms  involving  vertical  resonances  (e.g.,  Pfenniger  &  Friedli 
1991)  are  not  likely  to  be  effective,  since  the  single-orbit  res¬ 
onance  condition  is  likely  to  be  destroyed  by  the  presence  of 
collective  instabilities,  as  our  earlier  analysis  showed.  In  fact,  the 
so-called  buckling  instability  attributed  to  those  vertical  resonances 
acts  very  fast  (Pfenniger  &  Friedli  1991)  and  would  have  resulted 
in  two  distinct  classes  of  disk  galaxies  (one  before  the  buckling, 
with  no  bulge,  and  one  after  the  buckling,  with  a  bulge),  whereas 
the  observed  bulges  come  in  all  variations  of  bulge-to-disk  ratio, 
and  there  is  no  signature  of  the  above  bimodality  which  would  have 
resulted  if  buckling  were  the  underlying  mechanism  for  bulge  for¬ 
mation.  The  buckling  observed  in  previous  three-dimensional 
A-body  simulations  is  apparently  the  result  of  the  initial  conditions 
chosen  in  these  simulations,  which  results  in  the  rapid  evolution  of 
the  initial  configuration  to  reach  a  new  equilibrium  quasi- steady 
state. 

Note  that  past  A-body  simulations  have  had  difficulty  re¬ 
producing  the  large  accretion  rate  required  for  the  cosmological 
transformation  of  the  Hubble  types  of  galaxies  (e.g.,  Zhang  1998, 
1999).  This  is  now  understood  to  be  mainly  due  to  the  fact  that 
these  simulations  used  a  much  smaller  number  of  particles  than 
the  number  of  stars  present  in  physical  galaxies  (106  rather  than 
1011).  This  produced  an  increase  in  the  initial  growth  rate  of  bar 
modes,  as  well  as  in  the  collisional  relaxation  induced  by  the 
collective  instability,  which  together  demand  that  a  rigid  halo/ 
bulge  component  be  imposed  (especially  for  the  2D  A-body  simu¬ 
lations)  in  order  to  artificially  suppress  the  violently  growing  bar 


modes,  so  that  its  resulting  rapid  heating  would  not  suppress  the 
appearance  of  the  spiral  modes  (Ostriker  &  Peebles  1973).  How¬ 
ever,  the  suppression  of  the  density  modes  with  large  growth  rates 
also  means  that  in  these  simulations  we  have  only  simulated  the 
class  of  spiral  patterns  which  have  small  perturbation  amplitudes 
(the  density  perturbation  amplitudes  in  the  simulated  spirals  are 
usually  around  15%-20%,  and  the  potential  perturbation  am¬ 
plitudes  are  around  5%).  The  small  wave  amplitude  in  the  sim¬ 
ulated  spirals  is  the  main  reason  for  the  small  accretion  rate 
observed  in  A-body  simulations,  a  conclusion  already  reached  in 
Zhang  (2003).  These  facts  reflect  the  inherent  difficulty  in  using  a 
small  number  of  particles  to  simulate  a  highly  nonlinear  struc¬ 
ture  which  has  a  large  perturbation  amplitude.  Physical  galaxies 
have  significantly  more  particles  than  used  in  the  simulations  so 
far  and  thus  can  sustain  sufficiently  large  arm  /interarm  contrast 
for  their  density  waves  to  result  in  the  kind  of  secular  evolution 
rates  observed  for  cluster  and  field  galaxies. 

We  also  comment  here  that  in  the  past  other  researchers  have 
attempted  to  use  resonance  broadening  of  a  time-varying  density 
wave  pattern  to  reconcile  the  observed  secular  evolution  be¬ 
havior  of  particle  trajectories  found  in  A-body  simulations  with 
the  LBK72  conclusion  of  no-wave  basic-state  interaction  except 
at  resonances  (see,  e.g.,  Sellwood  &  Binney  2002  and  references 
therein).  We  point  out  that  the  secular  evolution  process  dis¬ 
cussed  in  this  work  is  not  a  direct  result  of  the  time-fluctuating 
nature  of  the  density  wave;  in  fact,  the  mass  accretion  rate  ex¬ 
pression  (eq.  [9])  becomes  exact  precisely  at  the  quasi-steady 
state  of  the  wave  mode  (Zhang  1998).  Indeed,  there  are  micro¬ 
scopic  time  fluctuations  of  the  overall  potential  in  any  A-particle 
disk,  but  such  microscopic  fluctuations  are  not  an  integral  part  of 
the  resonance-broadening  analysis.  (That  is,  one  could  use  dif¬ 
ferent  numbers  of  particles,  different  initial  conditions,  etc.,  in 
these  simulations  and  still  obtain  the  same  mode  and  secular 
evolution  rate,  as  long  as  the  basic  state  chosen  is  the  same.  The 
only  change  with  the  use  of  different  particle  numbers  is  the  re¬ 
sulting  heating  rate.  But  the  mass  accretion  rate  and  angular 
momentum  exchange  rate  remain  the  same.  See  Zhang  [1998], 
§  5.2,  for  further  details  of  these  tests.)  Another  indication  that 
the  secular  evolution  and  collective  dissipation  behavior  we  dis¬ 
cussed  are  unrelated  to  the  broadening  of  resonances  is  that,  as 
recently  found  by  Sellwood  &  Binney  (2002),  the  resonance¬ 
broadening  approach  produces  radial  migrations  of  particles  in 
an  overstable  disk  subject  to  transient  spiral  perturbations  that  are 
exactly  opposite  in  direction  to  what  is  predicted  in  secular  evo¬ 
lution  theory,  i.e.,  to  what  was  also  observed  in  A-body  simulations 
containing  unstable  global  modes.  The  resonance-broadening  ar¬ 
guments  thus  clearly  failed  to  capture  the  essential  physics  of  the 
secular  evolution  process  in  physical  galaxies  containing  global 
density  wave  modes. 

4.4.  Phase  Shifts,  Torques,  and  Angular  Momentum  Exchanges 

In  this  subsection  we  address  several  issues  related  to  the  phase 
shifts,  torques,  and  angular  momentum  transports  and  exchanges 
in  disk  galaxies,  which  underlie  the  approach  we  used  for  deter¬ 
mining  the  corotation  radii  and  calculating  the  mass  accretion 
rates. 

There  are  actually  two  different  meanings  when  we  use  the 
word  “torque”  in  the  context  of  spiral  waves.  There  is  the  torque 
expression,  which  (when  multiplied  by  2irr)  is  the  total  torque 
experienced  by  the  matter  within  a  unit- width  annulus  due  to  the 
torquing  of  the  wave  [we  called  this  T(r)  in  Zhang  (1996),  or 
—T\(r)  as  defined  in  Zhang  (1998);  here  the  negative  sign  in 
front  of  T\{r)  signifies  the  fact  that  T\{r)  represents  the  torquing 
in  the  —  ^-direction  by  the  wave  on  the  disk  matter  inside 
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corotation],  which  is  effectively  due  to  the  rest  of  the  disk  matter 
not  in  the  annulus.  There  is  also  the  torque  coupling  integral  that 
LBK72  used  in  their  paper  [which  can  be  further  broken  up  into 
the  gravitational  torque  coupling  integral  Cg{r)  and  advective 
torque  coupling  integral  Ca(r)],  which  expressed  the  torquing  of 
the  inner  disk  on  the  outer  disk,  or  vice  versa,  across  a  circular 
boundary  at  a  given  radius  r  of  the  disk.  The  radial  gradient  of 
LBK72’s  gravitational  torque  coupling  integral  dCg/dr  is  equal 
to  our  torque  T \  (r)  in  the  linear  regime  of  the  wave  mode  (here 
the  sign  convention  works  such  that  the  positive  dddr ,  which  sig¬ 
nifies  an  angular  momentum  drain  from  an  annulus,  is  related  to 
the  -^-direction  of  the  torque  of  the  wave  on  the  disk  matter).  At 
the  quasi- steady  state  of  the  wave  mode  it  is  actually  the  radial 
gradient  of  the  sum  of  LBK72’s  gravitational  and  advective  torque 
coupling  integrals  that  equals  the  torque  expression  T\  (r)  (Zhang 
1998, 1999).  Incidentally,  the  LBK72  total  torque  coupling  integral 
C(r)  =  Ca(r)  +  Cg(r)  is  the  same  as  the  total  outward  angular  mo¬ 
mentum  flux  across  a  given  radius  r,  if  we  are  dealing  with  a  spi¬ 
ral  density  wave. 

The  significance  of  these  relations  between  the  LBK72  torque 
coupling  integrals  and  our  torque  integral  Tfr)  is  that  LBK72’s 
total  torque  coupling  C(r)  needs  to  be  bell-shaped  in  order  for  the 
relevant  modes  to  be  able  to  spontaneously  and  homogeneously 
emerge  from  an  originally  featureless  disk,  and  this  bell-shaped 
torque  coupling  must  correspond  to  a  radial  distribution  of  T \  (r) 
[which  is  itself  proportional  to  sin  (m0o)>  where  0O  is  the  potential- 
density  phase  shift],  which  has  a  positive  hump  inside  corotation 
and  a  negative  hump  outside  corotation,  exactly  the  same  as  the 
distribution  of  the  phase  shift  we  have  determined  in  this  work. 
The  change  of  sign  of  the  phase  shift  at  corotation  can  be  shown 
to  be  a  direct  consequence  of  a  mode’s  capability  to  spontaneously 
emerge  out  of  an  initially  axisymmetric,  differentially  rotating 
disk.  This  is  because  in  a  differentially  rotating  disk,  the  density 
wave  mode  formed  has  negative  angular  momentum  density  (with 
respect  to  the  axisymmetric  state  of  the  disk)  inside  corotation  and 
positive  angular  momentum  density  outside  corotation.  In  order 
for  the  wave  mode  to  be  able  to  spontaneously  emerge,  angular 
momentum  has  to  be  taken  out  of  every  annulus  of  the  inner  disk 
inside  r(CR)  and  deposited  onto  every  annulus  of  the  outer  disk 
outside  r(CR)  (and,  unlike  what  LBK72  originally  thought,  this 
removal  and  deposition  must  happen  only  at  resonances).  This 
fact  also  explains  why  the  spontaneously  growing  type  of  mode 
had  to  be  of  spiral  shape  (even  the  bar  mode,  when  it  first  emerged 
from  an  A-body  disk,  was  invariably  shown  to  be  of  spiral  shape 
initially  and  then  gradually  evolved  into  a  straight  bar;  see,  e.g., 
the  simulation  of  Sparke  &  Sellwood  [1987]):  because  a  perfectly 
straight  bar  will  not  be  able  to  transport  angular  momentum  out¬ 
ward,  and  without  lowering  the  angular  momentum  density  in  the 
inner  disk  (up  to  the  corotation  radius  of  the  mode),  a  mode  will 
not  be  able  to  emerge. 

The  above  refers  mostly  to  the  linear  modal-growth  regime. 
Then,  at  the  quasi- steady  state  of  the  wave  mode,  since  by  def¬ 
inition  the  wave  amplitude  can  no  longer  grow,  yet  the  bell-shaped 
torque  coupling  remains  as  long  as  the  mode  still  contains  some 
degree  of  skewness,  and  thus  the  angular  momentum  is  being 
continuously  transported  outward  (as  well  as  being  taken  out  of 
every  annulus  of  the  inner  disk  and  deposited  onto  every  an¬ 
nulus  of  the  outer  disk),  the  out-transported  angular  momentum 
cannot  come  from  the  wave  itself,  but  has  to  come  from  the  basic 
state;  therefore,  basic-state  evolution  (i.e.,  matter  inflow  and  out¬ 
flow)  is  a  natural  outcome  of  the  combined  facts  of  the  quasi¬ 
stationary  nature  of  the  wave  mode  and  the  continued  outward 
angular  momentum  transport  by  the  wave  mode  at  its  quasi¬ 
steady  state. 


Such  a  wave/basic-state  interaction  at  the  quasi- steady  state 
of  the  wave  mode  is  not  due  to  the  interaction  of  the  mean  wave 
potential  with  a  single  orbit  (since  the  conservation  of  the  Jacobi 
integral  in  the  corotating  frame  of  the  mode  forbids  the  single  or¬ 
bit  to  lose  energy  and  angular  momentum),  but  has  to  be  through 
the  collective  instability  at  the  spiral  arms,  with  its  integral  ex¬ 
pression  given  by  equation  (1).  The  effect  of  the  torque  integral 
T (r)  on  a  group  of  stars  has  been  quantitatively  confirmed  through 
A-body  simulations  (Zhang  1998,  §  5.2.2;  see  especially  the  text 
in  relation  to  Fig.  10  in  that  paper). 

In  regard  to  the  “solid  body”  appearance  of  the  wave/basic- 
state  interaction  described  by  equation  (1),  despite  the  seeming 
randomness  of  each  star’s  orbital  motion  with  respect  to  the  wave, 
we  want  to  point  out  the  difference  in  torque  action  in  three  types 
of  scenarios:  (1)  A  true  solid  body,  as  in  the  case  of  a  curved 
sword  or  a  scimitar  (an  analogy  made  by  B.  Elmegreen  2006,  pri¬ 
vate  communication).  Here  we  have  the  gravitational  torque  and 
torque  coupling  (but  no  advective  torque  coupling,  since  the  par¬ 
ticle  configuration  is  not  in  motion),  yet  no  angular  momentum 
transport.  This  is  because  for  every  segment  of  the  sword  the 
torque  coupling  integrals  from  the  inner  and  outer  boundaries  of 
the  segment  exactly  cancel  each  other  (the  gravitational  torque 
coupling  still  has  a  gradient  as  in  the  case  of  a  spiral  mode,  but 
this  gradient  is  compensated  for  by  the  countering  torque  gradient 
due  to  the  molecular  bounding  forces).  That  we  can  have  torque 
and  countertorque  but  no  net  angular  momentum  exchange  in  this 
situation  is  similar  to  pushing  a  brick  wall:  one  can  apply  a  force 
and  experience  a  counterforce,  but  unless  the  wall  moves  one  will 
not  have  done  any  work  on  the  wall.  (2)  A  loose  association  of 
self-gravitating  particles  arranged  in  the  shape  of  a  curved  sword. 
In  this  case,  there  will  be  gravitational  torque  as  before  (and  there 
might  even  be  advective  torque  when  the  particles  start  moving), 
but  the  configuration  will  change  as  a  result,  and  there  is  indeed 
angular  momentum  transfer  between  the  inner  and  outer  portions 
of  this  loose  association  of  particles  (they  will  swing  in  opposite 
directions  to  straighten  out  the  curve  of  the  sword,  yet  the  total 
angular  momentum  of  the  system  remains  zero).  (3)  The  case  of 
spiral  modes  in  galaxies.  This  situation  is  unlike  that  of  either  (1) 
or  (2),  in  that  there  is  simultaneously  the  gravitation /advective 
torque  and  torque  coupling,  the  transfer  of  angular  momentum, 
and  the  maintenance  of  the  coherence  of  the  original  modal  con¬ 
figuration.  This  structural  stability  of  the  mode  in  the  face  of  per¬ 
turbing  influences  is  what  made  the  global  instability  patterns  so 
unique.  The  structural  stability  is  achieved  as  a  dynamical  balance 
between  the  spontaneous  growth  tendency  of  the  unstable  mode 
and  the  dissipative  damping  tendency  of  the  collective  viscous 
forces.  The  wave  acts  as  if  it  is  a  solid  entity  which  can  receive  and 
give  out  energy  and  angular  momentum,  yet  it  simultaneously 
possesses  a  soft  boundary  in  the  sense  that  the  constituting  stars 
which  support  the  wave  motion  swarm  in  and  out  of  the  wave 
crest. 

4.5.  The  Physical  Basis  Underlying  the  Validity 
and  Accuracy  of  the  Phase-Shift  Method 

The  validity  of  the  phase-shift  method  for  corotation  determina¬ 
tion  is  rooted  in  the  global  self-consistency  requirement  of  the 
wave  mode  (i.e.,  both  the  Poisson  equation  and  the  equations  of 
motion  need  to  be  satisfied  at  the  same  time).  Therefore,  the  ap¬ 
pearance  that  the  corotation  can  be  determined  from  the  Pois¬ 
son  equation  alone,  without  knowledge  of  the  kinematic  state 
of  the  disk  matter,  is  only  an  illusion.  In  actuality  the  global  self- 
consistency  between  the  Poisson  equation  and  the  equations  of  mo¬ 
tion  has  been  enforced  by  nature  itself  when  the  galactic  resonant 
cavity  filters  and  selects  the  set  of  modes  it  can  support.  The  phase 
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shift  resulting  from  the  Poisson  equation  is  determined  not  only  by 
the  spiral  pitch  angle  but  by  the  radial  density  variations  of  the 
nonaxisymmetric  density  wave  features  (Zhang  1996).  For  a  self- 
sustained  global  spiral  mode,  the  radial  density  variation  of  the 
modal  perturbation  density,  as  well  as  the  pitch  angle  variation, 
together  will  be  such  that  the  Poisson  equation  will  lead  to  the  zero 
crossing  of  the  phase-shift  curve  being  exactly  at  the  corotation 
radius  of  the  mode,  since  that  is  where  the  peak  of  the  bell-shaped 
angular  momentum  flux  is,  and  the  radial  gradient  of  the  angular 
momentum  flux  is  proportional  to  the  sine  of  the  phase  shift 
(Zhang  1998;  see  also  eq.  [C8]  of  Lin  &  Lau  1979).  Since  the  bell¬ 
shaped  angular  momentum  flux  distribution  is  what  allowed  the 
spontaneous  growth  of  the  mode  in  the  linear  regime,  only  modes 
with  density  distributions  which  result  in  the  correct  phase-shift 
distributions  are  naturally  selected  by  the  galactic  resonant  cavity. 

The  self-consistency  requirement  of  a  spontaneously  formed 
density  wave  mode  explains  not  only  why  we  can  use  the  Poisson 
equation  alone  to  determine  the  corotation  radius,  but  also  why  we 
can  use  the  continuity  equation  alone  to  determine  the  pattern 
speed  as  in  the  TW84  approach:  each  is  only  seeing  one  side  of 
the  coin,  but  since  both  sides  join  at  a  unique  boundary,  we  can  de¬ 
termine  the  circumference  of  the  coin  from  measurements  con¬ 
ducted  on  either  side. 

The  phase-shift  distribution  calculated  through  the  Poisson 
equation  would  be  the  same  whether  the  pattern  were  rotating  or 
not,  as  long  as  we  used  the  same  surface  density  distribution.  The 
phase- shift  method  is  thus  inherently  not  a  rotation  detection 
method.  Rather,  it  is  a  diagnostic  tool  for  characterizing  the  fea¬ 
tures  of  self-organized  modes  in  galactic  disks.  We  can  imagine, 
for  the  sake  of  argument,  that  a  mode  had  initially  reached  a 
quasi- steady  state  with  a  given  pattern  speed.  If,  suddenly,  we 
apply  a  forced  potential  of  exactly  the  same  shape  as  the  original 
modal  potential,  but  with  a  somewhat  faster  pattern  speed  than  the 
original  natural  pattern  speed,  the  pattern  density  will  tend  to 
follow  this  applied  potential’s  forcing  with  only  minor  changes 
in  shape.  The  phase-shift  method  will  not  be  able  to  predict  that 
the  rotation  speed  of  the  pattern  has  changed  in  this  case  or, 
equivalently,  that  the  corotation  radius  has  changed.  But  once  the 
forcing  is  removed,  the  pattern  will  quickly  recover  its  original 
pattern  speed.  The  reason  is  that  the  faster  pattern  speed  has  moved 
the  corotation  radius  from  its  original,  natural  position,  which 
previously  allowed  the  sign  of  the  angular  momentum  deposition 
by  the  wave  everywhere  across  the  disk  to  match  the  sign  of  the 
angular  momentum  density  of  the  wave  mode  itself.  This  match¬ 
ing  condition  is  violated  as  in  the  forced  case;  once  we  remove 
the  forcing,  the  faster  mode  will  be  damped  out  because  of  the 
wrong  sign  and  wrong  amounts  of  angular  momentum  being 
deposited  by  the  wave  mode  over  much  of  its  extent.  The  original 
mode  will  reappear,  since  it  is  naturally  amplified  in  the  galactic 
resonant  cavity  determined  by  the  basic  state.  Or,  to  put  it  in 
another  way,  the  spiral  shock  and  the  associated  collective  dissi¬ 
pation  tend  to  destroy  or  damp  the  wave  at  every  instant  of  time.  It 
is  only  through  the  continued  amplification  tendency  of  the  orig¬ 
inal  mode  in  the  galactic  resonant  cavity  that  a  quasi- steady  state 
can  be  maintained.  The  modes  that  satisfy  the  self-sustainability 
condition  have  a  special  appearance  that  matches  their  kinematic 
characteristics,  so  that  the  growth  tendency  balances  the  dissi¬ 
pation  tendency  to  reach  a  dynamical  equilibrium  at  the  quasi¬ 
steady  state. 

We  thus  see  that  as  a  general  method  for  corotation  determi¬ 
nation  the  visual  matching  of  the  forced  pattern  response  with 
that  of  the  observed  galaxy  morphology  is  not  likely  to  be  reli¬ 
able,  since  the  response  pattern  tends  to  look  almost  identical  for 
quite  a  wide  range  of  pattern  speed.  The  limitations  of  the 


passive  simulation  approach  become  all  the  more  prominent 
when  we  realize  that  realistic  galaxies  contain  nested  density 
wave  patterns  which  all  rotate  at  different  pattern  speeds,  as  we 
have  found  in  this  work. 

On  the  other  hand,  the  phase-shift  method  itself  is  only  ex¬ 
pected  to  be  reliable  for  cases  of  spontaneously  formed  density 
wave  patterns  that  have  reached  a  quasi-steady  state.  We  have 
found  empirically  from  analyzing  the  OSUBGS  that  the  early- 
type  disk  galaxies  usually  give  the  best  fit  between  the  predicted 
corotations  and  the  resonance  features,  likely  related  to  the  fact 
that  early-type  galaxies  have  had  the  longest  time  to  evolve  and 
reach  a  quasi- steady  state.  The  latest  Hubble  types  usually  have 
more  flocculent  patterns,  and  the  phase-shift  result  correspond¬ 
ingly  does  not  show  a  well-organized  pattern,  partly  reflecting 
the  fact  that  the  wave-mode  pattern  speeds  and  corotation  radii  in 
these  galaxies  are  themselves  poorly  defined.  The  intermediate 
Hubble  types  have  strong  and  open  spiral  arms,  but  the  results  of 
the  phase-shift  calculation  tend  to  be  more  sensitive  to  the  uncer¬ 
tainty  in  inclination  angle  deprojections  which  can  result  in 
inaccurate  spiral  pitch  angle  representations,  an  issue  we  will  ad¬ 
dress  further  in  the  sequel  paper.  The  intrinsic  accuracy  of  the 
phase-shift  method  should  be  high  if  we  are  dealing  with  an  ideal, 
spontaneously  formed  mode  that  has  reached  a  quasi- steady  state, 
and  if  we  have  accurate  knowledge  of  the  orientation  parameters 
of  the  galaxy.  In  practice,  however,  the  mode  may  not  be  in  a 
quasi- steady  state,  and/or  the  orientation  parameters  may  be 
uncertain.  Since  the  errors  are  likely  to  be  dominated  by  unknown 
systematics  rather  than  by  random  measurement  errors,  we  have 
not  assigned  error  bars  to  our  results.  Other  sensitivities/pitfalls  of 
the  method  will  be  discovered  as  it  is  applied  to  a  greater  number 
of  real  and  simulated  galaxies. 

We  now  comment  briefly  on  the  pattern  speeds  and  corotation 
radii  determined  from  previous  modal  calculations  and  their 
relation  to  the  results  of  our  current  study.  In  the  1 970s  and  1980s 
Lin  and  collaborators  constructed  self-consistent  spiral  modes 
both  through  solving  the  exact  fluid  set  of  equations  in  the  linear 
regime  and  through  the  asymptotic  approach  of  higher  order 
local  calculations  which  obtained  piecewise  continuous  wave 
trains,  and  then  joined  them  together  to  derive  the  quantum  con¬ 
ditions  for  the  pattern  speed  determination  (see,  e.g.,  Lin  &  Lau 
1979  and  references  therein;  Bertin  et  al.  1989a,  1989b  and  ref¬ 
erences  therein).  These  two  approaches  were  found  to  agree  with 
one  another  as  long  as  one  deals  with  slow-growing  modes.  In 
both  approaches,  the  quantum  condition  for  the  mode  was  derived 
by  enforcing  the  global  self-consistency  requirement  between  the 
Poisson  equations  and  the  equations  of  motion  (or  the  so-called 
Poisson-Euler  set  in  the  fluid  approach).  In  the  asymptotic  ap¬ 
proach,  this  global  self-consistency  was  acquired  in  a  piecewise, 
locally  self-consistent  manner  and  by  throwing  away  the  so-called 
out-of-phase  terms  (which  exactly  corresponds  to  the  potential- 
density  phase  shift  we  discuss  in  this  work).  The  fact  that  this  latter 
approach  leads  to  results  similar  to  the  “exact”  numerical  solu¬ 
tions  which  accurately  contain  all  the  phase-shift  information  is 
because  the  asymptotic  approach  only  solved  the  weakly  growing 
kind  of  modes  (i.e.,  either  very  small  pitch  angles,  as  in  tightly 
wound  early  Hubble  type  galaxies,  or  very  large  pitch  angles,  as 
in  very  straight  bars,  again  appearing  in  early  Hubble  type  gal¬ 
axies;  both  of  these  cases  give  small  phase  shifts,  since  the  phase 
shift  is  largest  for  pitch  angles  close  to  45°  and  is  zero  for  pitch 
angles  of  0°  and  90°).  Therefore,  the  asymptotic  approach  worked 
precisely  for  situations  where  the  phase  shift  was  small,  so  the 
resulting  small  out-of-phase  terms  (which  reveal  the  effect  of 
potential-density  phase  shift)  can  be  thrown  away.  For  the  fast¬ 
growing  “violently  unstable”  modes,  which  have  a  more  open 
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spiral  morphology  (which  resembles  intermediate  Hubble  type 
galaxies),  the  exact  numerical  solutions  of  linear  global  modes 
do  contain  the  phase-shift  distribution  that  changes  sign  at  coro¬ 
tation  (C.  C.  Lin  1995,  private  communication),  and  it  is  precisely 
this  phase-shift  distribution  in  the  open  spiral  modes  which 
leads  to  the  large  growth  rate  of  the  mode  in  the  linear  regime 
and  to  the  rapid  azimuthal  steepening  of  the  sinusoidal  wave 
profile  into  narrow  shock  fronts  in  the  nonlinear  regime  (Zhang 
1998). 

Realistic  galaxies  also  often  contain  companions  and/or  have 
gone  through  major/minor  mergers.  In  future  studies  we  will 
explore  the  role  of  environment  on  the  reliability  and  detailed 
distribution  of  phase  shifts  in  galaxies.  We  comment  here  that  if 
the  galaxy  disks  involved  are  unstable  to  the  formation  of  intrinsic 
modes,  the  effect  of  an  interaction  often  only  serves  to  excite  a 
large-amplitude  mode  which  has  the  same  properties  as  the  in¬ 
trinsic  mode  (Zhang  et  al.  1993);  i.e.,  the  intrinsic  modes  in  un¬ 
stable  disks  can  be  rather  robust  during  modest  galaxy  interactions. 
Minor  mergers  potentially  can  result  in  a  change  of  the  basic 
state  and  thus  the  kind  of  unstable  modes  admitted.  The  phase- 
shift  method  should  be  applicable  in  this  case  as  soon  as  the  entire 
configuration  has  settled  into  the  new  quasi- steady  state.  Major 
mergers,  on  the  other  hand,  often  result  in  damages  that  are  so 
severe  that  the  galaxy  morphology  is  no  longer  described  by  grand- 
design  density  wave  patterns,  so  the  phase-shift  method  will  not 
apply.  For  certain  fortuitous  cases  such  as  NGC  4622  the  method 
could  still  work  even  for  transient  waves,  but  in  general,  its  validity 
has  to  be  checked  on  an  individual  basis  for  any  galaxy  involved  in 
interactions. 

5.  CONCLUSIONS 

In  this  work  we  have  developed  and  verified  a  general  approach 
for  determining  the  corotation  radii  of  density  waves  in  disk  gal¬ 
axies,  which  applies  to  any  disk  galaxies  whose  density  wave 
modes  have  reached  a  quasi- steady  state,  a  condition  empiri¬ 
cally  found  to  be  the  case  for  the  majority  of  the  nearby  disk 
galaxies. 


Compared  to  other  approaches  developed  so  far,  the  new  ap¬ 
proach  for  corotation  determination  has  a  number  of  advantages: 
(1)  It  is  relatively  insensitive  to  star  formation,  MIL  variations, 
and  vertical  scale  height  assumptions,  at  least  at  the  level  these 
problems  are  present  in  the  NIR.  (2)  It  can  be  applied  to  face-on 
galaxies;  the  orientation  of  the  bar  has  little  impact.  (3)  It  can  be 
effectively  applied  to  all  Hubble  types,  at  least  those  with  a  disk 
shape.  (4)  The  radial  extents  of  the  nested  resonances  and  mul¬ 
tiple  pattern  speeds  can  be  determined  in  a  model-independent 
fashion,  and  in  many  cases  the  central  resonances  can  be  deter¬ 
mined  with  “superresolution”  due  to  the  nonlocal  nature  of  the 
potential  calculation.  (5)  The  application  of  the  method  can  use 
existing  databases  of  NIR  galaxy  images  without  the  need  for 
significant  investments  in  new  telescope  time.  Furthermore,  the 
simplicity  of  the  application  of  the  method  allows  the  statistical 
analysis  of  the  correlation  of  the  calculated  corotation  radii  with 
the  observed  morphological  features  of  galaxies.  (6)  The  con¬ 
nection  between  the  potential-density  phase  shifts  and  the  sec¬ 
ular  evolution  rates  of  galaxies  facilitates  the  study  of  the  galaxy 
morphological  and  kinematic  features  in  the  context  of  the  cosmic 
evolution  of  galaxies.  This  provides  another  way  of  quantifying 
the  impact  of  bars  and  spirals  on  disk  galaxies. 
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